Under Pressure: Star Clusters and the Neutral Hydrogen 

Medium of Tidal Tails 
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ABSTRACT 

Using archival data from ATCA, WSRT, and the VLA, we have analyzed the Hi emission 
of 22 tidal tail regions of the Mullan et al. sample of pairwise interacting galaxies. We have 
measured the column densities, line-of-sight velocity dispersions, and kinetic energy densities on 
~kpc scales. We also constructed a tracer of the line-of-sight velocity gradient over '-^lO kpc 
scales. We compared the distributions of these properties between regions that do and do not 
contain massive star cluster candidates (My < -8.5; ~ 10^-10^ Mq as observed in HST WFPC2 
VI data). In agreement with Maybhate et al., we find that a local, ~kpc-scale column density of 
log Nhi ^ 20.6 cm~^ is frequently required for detecting clustered star formation. This Hi gas 
also tends to be turbulent, with line-of-sight velocity dispersions ctios ~ 10-75 km s""'^, implying 
high kinetic energy densities (log Ske > 46 erg pc~^). Thus high Hi densities and pressures, 
partly determined by the tail dynamical age and other interaction characteristics, are connected 
to large-scale cluster formation in tidal tails overall. Lastly, we find that the high mechanical 
energy densities of the gas are likely not generally due to feedback from star formation. Rather, 
these properties are more likely to be a cause of star formation than a result. 

Subject headings: galaxies: ISM, galaxies: interactions, galaxies: star clusters, radio lines: galaxies 
1. Introduction 
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Tidal tails are common products of disk galax y 
interactions (e.g., lToomrelll977l : ISchweizeiill978l) . 
These extended, pronounced structures begin 
their lives shortly after the initial encounter of 
their progenitor galaxies, drawn from a com- 
bination of ISM and stars of these systems' 
outer regions. While relatively short-lived (<1 
Gyr), tidal debris exhibits an extraordinary va- 
riety of morphological and kinematic phenom- 
ena. These range from long, thin tails, to dif- 
fuse plumes, to kinematically distinct tidal dwarf 
galaxies (TDGs), all characteristic of the numer- 
ous interactions that populate the Local Volume 
(iDuc et al.ll200l . 

Tidal tails are composed of ~10^ Mq of 
gas — and some small underlying population of 
stars — pulled from the outer regions of their host 
galaxies by gravitational torques and the tidal 



field indu ced by an encoun ter with another sys- 
tem (see iBournaudI I2OIII for a reyiew). Al- 
though there are wide yariations in metallic- 
ity by galaxy, studies of non-interacting disk 
galaxies indicate that gas from these regions 
prior to the inte raction are genera lly Hl-rich and 
metal/H2-poor (JBigiel et al.l 120101 ). The neutral 
hydrogen medium is thus the de facto tracer 
of the kpc-scale ISM in tidal tails, and has a 
long history of obseryati on 



of molecular g a s ayailable for collapse (e.g. 



19791: iSimkinet al 



[Ml 



yan der Hulst 



19871: lApp leton et al.' 'l98 



Yun et all Il994t iHorellou fc Koribalski 20071 
Sengupta et al.ll2012[ ). Decades of study have re- 
vealed that the Hi in tidal debris is a generally 
inhomogeneous, kinematically disturbed medium, 
subject to a variable tidal potential throughout 
the course of the interaction that can be both 
locally "extensive" (pulling materi al apart) and 
"compressive" (pushing it together; iRenaud et al. 
I2OO9I) . 

The importance of Hi column density, Nhi, 
in star formation has observational precedent 
in many extragalactic environments. In a self- 
gravitating gas disk, stability analysis s hows there 
is like ly a threshold Nhi that permits it (JKennicutt 
Il989l ). This may be true across a variety of cosmic 
venues. On ~500 pc^ scales within dwarf irregu- 
lar galaxies, this threshold may be ftilO ^^ cm~^ 
(|Skillman fc Bothunl[l98i ISkillmanlll987l ). More- 
over, UV emission — a key diagnostic for recent 
star formation in the past ^100 Myr — is observed 
where the Hi colu mn dens ity reaches or e xceeds 
10^0 cm-^ (iNeff et al. 2005; Thilkcr et al 



i2.5 X 



2005t iGil de Paz et al.ll2005l) . This is consistent 



with ISchavd (|2004l) , who determine that 



a mini- 
mum surface density of ss3 Mq pc~^ is necessary 
to first form a cold neutral phase, with some dis- 
persion from metallicity and the interstellar radi- 
ation field. The production of the cold molecu- 
lar hydrogen necessary for gravitational collapse 
and star formation would then occur, the amount 
of which is controlled largely by the Hi density 
iKrumholz et al.ll2009t ). The ratio of H2 to Hi sur- 
face densities increases almost linearly with inter- 
stellar pressure and the Hi densities that promote 
it, and saturates at a mas s density of ^10 M, 



££! 



2004; 



1 
Lei 



Wong & Blitz' '2002'; 'Blitz 
Blitz & Rosolowsky 2006; 
Star 



20081: iLerov et al.l I2OO8I) 

then scale with respect to the density and amount 



Rosolowsk 
; Walter ct al 
formation may 



i 



Liu et al. 2011 and references 



iKennicutd Il998 : 
therein). 

In tidal tails, large distances (<100 Mpc) 
and low surface brightnesses can make observing 
this star formation difficult. Compact star clus- 
ters and their populations are therefore conve- 
nient surrogates for this task; they can be con- 
sidered high-signal, luminous tracers of the lo- 
cal star formation history of their en vironments 
(|Konstantopoulosll2010t Ide Griidbofoh . The most 
massive ('^lO^-lO^ -^o) of these star clusters or 
cluster can didates have been observed in some 
tidal tails (Ide Griis et all l2002l: iTran et al.ll2003t 



201 ll : iKnierman et al.l l2003l ) . This implies that 



Bastian et al.l l2005l). but n ot all (JMuUan et al 



the conditions required for extensive star cluster 
formation (and the detection of the most mas- 
sive clusters by the size-of-sample effect; e.g., 
Gieles fc Bastianll2008[) are not always present in 



these environments. Consequently, tails may oc- 
cupy the edge of the physical parameter space that 
allows these structures to form and/or survive. 

The content and character of the observable Hi 
medium may have a role to play in star cluster for- 



Mavbhate et al.l (|2007l ) examined a sam- 



mation. 

pie of tidal debris from late-stage (>400 Myr old) 
interactions and determined that a local, ^kpc- 
scale Hi column density of wlO^"'^ cm~^ or higher 
is required for finding significant populations of lu- 
m inous star clu s ters. T his is similar to the results 
of I Walter et al.l ( 20061) . who find that Hii regions 



are only observed where log Nhi > 21.0 cm on 
200-pc2 scales for the tidal arms of NGC 3077. 



However. iMavbhate et al.l (|2007l ) stipulate that 
Nhi is not a sufficient condition for finding bright, 
massive star clusters. They cite the debris of 
NGC 3921 as an environment where a non- 
negligible fraction of neutral hydrogen exists above 
their fiducial threshold value, but massive star 
clusters are unlikely to exist. Other proper- 
ties related to the kinematics of the Hi medium 
may also be important. For instance, disrup- 
tive shear might influence a star-forming environ- 
ment to "prefer" forming asso ciations over bound 



clusters (|Weidner et al.ll2010l) . This may be ob- 



served in__tlie_ceiitral region s of the Antennae, 
where IWhitmore et al.l ( 1999f) do not find massive 
clusters where the line-of-sight velocity gradient is 
high. 



Morphological asymmetries in the optical and 
Hi debris components also indicate that other, 
large-scale physics may influence the Hi-star clus- 
ter connection. Dust obscuration, ram-pressure 
stripping by the galaxy interaction, and photoion- 
ization from a central galactic superwind or locally 
by a previous generation of efficient star forma- 
tion ( Hibbard et al.l 120001 ) may additionally con- 
fuse the issue. Some of these phenomena, along 
with the bulk properties of the tidal tail, will be 
dependent on age, initial galactic orientations and 
speeds, and other propertie s of the interactions 



themselves (JDuc et al.ll2004[ ) 



The inconsistent production of star cluster pop- 
ulations in tidal tails may also be partially due 
to the kinematic character of the ^kpc-scale Hi 
medium. Under "normal" circumstances in many 
parts of isolated galaxies, the Hi medium generally 
consists of a mixture of cold (T ~ 100 K) Hi clouds 
embedded in a warm Hi medium (T ~ 8000 K) 
with characteristic thermal l i newidths ^6-8 km 
s-^ (iMcKee fc Ostrikeilll977t IWolfire et al.lll995l: 
IWolfire et al.l 120031 ) . Past the canonical optical 
boundary i?25, the warm ISM may be heated by 
UV radiation from th e interior gala xy or the extra- 
galactic background (ISchavd 120041) . In tidal tails, 
however, observed linewidths often exceed ther- 
mal values. These motions are typically attributed 
to turbulence in the gas, which may be gene rated 



by fe edback from star formati on (Tamburro et al 



2009 ) or magnetic effects fe.g.. lSellwood fc Balbus 



19991 ). Perhaps most importantly, turbulence can 



also be produced by the large-scale shocks and 
gravita tional instabiliti es triggered by the inter- 
action (JBournaudlbOlOl ). 

The pressures exerted by turbulent warm Hi 
and the reservoir of high-density gas they pro- 
duce may therefore engender the conditions nec- 
essary for pockets of the cold neutral medium 
(jWolfire et al.1120031 ) to exist on unobservable ^pc 
scales. If provided ample dust, these regions 
become sites for H2 production and subsequent 
star formation. Observationally, interacting sys- 
tems do indeed appear to generally contain more 
molecular gas than t heir quiescent counterparts 



( Casasola et al.l 120041 ) . If gas becomes too tur- 



bulent, however, it may interfere w ith its own 
cooling processes (jWolfire et al.ll2003l) and prevent 
high-density molecular gas from either forming, 
or cooling efficiently if already present. For ex- 



ample, iGuillard et al.l (|2012l) detect highly turbu- 
lent CO (with velocity dispersions >100 km s^^) 
in coUisionally shocked regions of Stephan's Quin- 
tet. These regions also have low PAH-to-CO sur- 
face luminosities, indicating low star formation ef- 
ficiency from an evidently high gas heating rate. 
Moreover, sufficient Hi density must be required 
to build enough H2 to self-shield from photodisas- 
sociation by the extragalactic background or rel- 
atively nearby stellar sources. Thus, it is likely 
that star formation is contingent on Hi density as 
well as kinematics, dictated by large (~kpc)-scale 
phenomena. 

Interstellar pressure not only enhances star 
formation in general, but it may also help 
dictate the types of stellar "packaging" that 
are produced. In tidal tails, the presence 
of dense, turbulent gas may conceivably en- 
courage the form a tion of tidal dwarf ga.laxies 



1993 : 



(|Elmegreen 
or other bound 
clusters. High pressures 
ter from disr uption in its 



Elmegreen fc Efremovl Il997t ) 

structures like compact star 

help prevent a clus- 

early gas expulsion 



phase (e.g.. IWhitmore. Chandar. fc Fal 



muision 

J Hool 



Bastian fc Gielesll2008t and references therein). If 
star formation depends on gas density, then high 
pressures promote clustered star formation at a 
rate that increas es as the local star forma tion rate 
(SFR) increases ("Larsen fc Ri chtleiJlJOOol ). In ad- 
dition, [Elmegreen (2008) finds that bound clus- 
ters can form at lower average densities at high 
Mach numbers than at low Mach numbers; thus 
the threshold density needed to form clusters is 
effectively reduced in highly turbulent media and 
clusters can become a favored "mode" of star for- 
mation. 

Tidal tails appear to offer a unique environment 
where the individual Hi- related variables that may 
influence star cluster formation can be systemati- 
cally tested. In this paper, we explore the complex 
relationship between galaxy interactions, their ef- 
fects on the Hi medium in their tidal tails, and 
the star clusters that are selectively found within. 
This study is the first detailed investigation of the 
kpc-scale column densities and kinematic charac- 
ter of the neutral hydrogen medium and its rela- 
tionship to star cluster populations in this type of 
extragalactic environ ment. We extend the sa mple 
of merging galax i es of [ Mavbhate et al.l (|2007l ) and 



esofE 
2003f l 



Knierman et al.l ()2003l ) with archival Very Large 



Array ( VL A) , Australia Telescope Compact Array 
(ATCA), and Westerbork Synthesis Radio Tele- 
scope (WSRT) 21-cm observations of the tidal tail 
sample studied by iMuUan et alj (|201l[ ). The 22 
tails total span a wide range of dynamical ages, 
merging mass ratios, and previously established 
star cluster population sizes and characteristics. 
An overview of the Hi observations and reduc- 
tion is provided in Section [21 along with selection 
of star cluster candidate positions and photome- 
try from optical data. In Section [3l we present 
our methods of measuring various Hi quantities, 
e.g., column density and velocity dispersion. We 
examine the local Hi characteristics of the tails 
with respect to the cluster populations in Section 
m We present our results in Section [5l and discuss 
trends and differences in the H i content among es- 
tablished types of interactions. In this section we 
also examine turbulence and its role as a cause or 
result of star formation. Lastly, we offer conclu- 
sions and final remarks in Section [HI 

2. Observations 

2.1. The Tidal Tail Sample and Their Star 
Cluster Candidates 



We obtained the HST WFP C2 data for 
the tidal tail regions s tudied bv iMuUan et a.1. 
hereafter Mil, Mavbhate et al.l 



(20111 



and iKnierman et al.l ( 20031) . 



(120071) . 
These include 
contributions from cycles 7 (GO-7466) and 16 
(GO-11134), with imaging performed in the 
F555W/F814W and F606W/F814W bands, re- 
spectively. The systems that produced these tails 
were originally selected to represent a large pa- 
rameter space of observable and dynamical char- 
acteristics of galaxy interactions present in the lo- 
cal universe, from optical/Hi properties, to pro- 
genitor mass ratios, to interaction ages (Mil). 
In this study, we use the same naming conven- 
tion for debris regions as these publications; e.g., 
NGC 1487E and NGC 1487W refers to the east- 
ern and western tails of NGC 1487, respectively, 
while the three WFPC2 pointings of NGC 4038/9 
are denoted by NGC 4038A, NGC 4038B, and 
NGC 4038C. 

We also acquired the source catalogs of J2000 
positions and photometry of co mpact objects de- 
tected in lMavbhate et al.l (|2007l ) and Mil; the lat- 
ter contains sources detected by iKnierman et al 



( 20031 ). Source photometry had been transformed 
from F555W/F606W and F814 W to standard 



V an d / bandpasses with the iHoltzman et al 



(|l995l ) prescriptiorl^ in these studies to facil 
itate comparison with the earlier photometry 
published in standard magnitudes. Photometry 
had a lso been corrected for Galact ic extin ction 
us ing Schleeel. Finkbeiner. &: Da.yisl ( 19981) and 



Girardi et al. (120081): see Mil and Knierman et al. 



(|2003[ ) for details. Objects detected and measured 
with this methodology were defined as star clus- 
ter candidates (SCCs) if they met the following 
color-magnitude criteria: 



• My < -8.5. Many studies use a sim- 
ilar nia£nitude__cutoff__(Mi^^^^8 to ■ 
9: IKnierman et al.l [20031: iKniermanI jJOO 



Schweizer et al.lll99alWhitmore et al.lll99 
because of detection limits and the threat 
of contamination from non-cluster sources at 
fainter magnitudes. Source detection within 
the tail regions of the Mil sample was previ- 
ously found to be more than 50% complete 
at My < -8.5, so we implement this crite- 
rionhere. Mil also reaffi rmed the findings of 
Efremov &: Ivanovl ( 19871 ): i.e. that this mag- 



nitude limit sufficiently protects against con- 
tamination from single main sequence and 
post main sequence stars. 

V — I < 1.0. This allows for a wide range of 
possible metallicities, extinctions, and ages 



in ou r cluster population 



1993; 



(e.g. . 

6 Tonrv 



Aihar. Blakeslee. 
Whitmore et al. 1995; Kundu & Whitmor 
20011 ). This V — I limit was also selected to 



Peterson 
'199 



reflect uncertainties o r additional reddening 
from dynamical (IPortegies Zwart et al 



2010l ). stochas tic (iMaiz Apellanizj l2009l). 
and disruptive (JAnders et al.ll2009l ) effects. 



For similar reasons, Mil and iMavbhate et al 
(|2007l ) used V — I < 2.0 as a color limit in their 
work. In practice, comparing the completeness 
curves of in-tail vs. out-of-tail areas reveal rel- 
atively limited tail extinctions (Ay ~ 0.5; Mil), 
and stochastic effects should be restricted toV — I 



^Mll employed updated methods for recent data 
and instrumental characteristics documented in 
|http : //purcell ■ as ■ arizona ■ edu/wf pc2_calib/2608_07_19 . html | 



changes at the <5% level (JAnders et alj l2009t ) 
for the Mv < -8.5, -10^-10^ Mq clusters we 
study. Furthermore, an examination of the Mil 
sources and their color-magnitude diagrams re- 
veals that a more stringent cutoff helps reduce the 
number of out-of-tail and presumably background 
sources detected, which tend to be redder {V ~ I 
> 1). These results are substantiated by the mul- 
tiwavelength observa tions and selection criteria of 
Tranchoet aLl(|2012[ ). 



2.2. Hi Calibration and Reduction 

Table 1 lists the tails from 



Mil, 



Knierman et al.l (|2003f ). and iMavbhate et al.l 
(|2007t) for which archival Hi data exist (22 
tails total). We include the dynamical ages 
and mass ratios of these interacting systems for 
later reference. Our mass ratio values and tail 
ages are taken directly from Mil and sources 
therein. The mass ratio is defined as M1/M2, 
where Mi and M2 are the respective masses of 
the perturbing and main galaxies. Upon further 
review, we inverted the reported mass ratios of 
NGC 2782E and NGC 4747 (to 4 and 20) to 
reflect the understanding in the literature that the 
tail regions examined here are likely pulled from 
the smaller-massed galaxies of their interactions. 
Neutral hydrogen data for NGC 1614, 
NGC 2782, NGC 2992/3, NGC 2444, and 
NGC 2535 were taken from the VLA archive^ 
The data were reduced in AIPS using the stan- 
dard procedures detailed in the AIPS CookbooH^- 
Here we provide a brief outline of this method- 
ology. In addition to the program source, a 
bright primary (gain) calibrator and a secondary 
(phase) calibrator were observed. The primary 
calibrator was usually observed at least once in 
the observing program and the secondary calibra- 
tor observed every approximately 40 minutes to 
an hour of time on source. Using the channel-0 
data, the flux density of the primary calibrator 
was found using setjy. Then the amplitudes and 
phases of both the primary calibrator and the 
secondary calibrator were determined using calib. 
The flux density of the secondary calibrator was 



determined from the primary calibrator using the 
task getjy. The calibration was applied to the 
source using deal. The channel-zero calibration 
solutions were copied to the spectral line data. 
Finally, the spectral line data was bandpass 
calibrated in bpass using the primary calibrator. 

The calibrated data were test imaged in imagr 
and the line-free channels determined. Using this 
information, the data were continuum subtracted 
using uvlsf. The final, continuum subtracted im- 
ages were produced in imagr, the data were im- 
aged with var^dn^j^lliSSofthe "robust" weighting 
parameter ( Briggs et al.lll999l ). with values from 
0-5. A larger robust parameter enhances detec- 
tion of diffuse Hi emission at the expense of angu- 
lar resolution. In most cases, the diffuse emission 
was resolved out of the high resolution image. In 
the case of two overlapping intermediate frequen- 
cies (IFs), the images in the overlapping channels 
were averaged and a combined cube formed. 

We acquired a reduced data cube for NGC 4747 
from the Westerbork Synthesis Radio Telescope 
(WSRT) web archivcQ. These data were origi- 
nally taken, calibrate d, and reduced as a part o f 
the WHISP program (jvan der Hulst et al.lBoOll) . 
Here, we use the cube generated to «60" reso- 
lution, which we found was an appropriate bal- 
ance between spatial resolution and detecting low- 
signal Hi in the tidal regions of NGC 4747. 

The neutral hydrogen data for NGC 6872 and 
NGC 1487 were obtained from the ATCA archives. 
The NGC 1487 da ta were processed using Miriad 
(jSault et al.lll995l ) for all steps, including the cal- 
ibration (with the tasks uvflag, blflag, and mfcal) 
and removal of continuum (with the task uvlin). 
The primary calibrators were 1934-638 and 1151- 
348, while the secondary was 0332-403. We pro- 
duced images by combining the array configura- 
tions listed in Table 1 — using standard tasks in- 
vert, clean, and restor. The data for NGC 6872 
were calibrated with a combination of AIPS and 
Miriad. AIPS was used to flag (spflag) and im- 
age the data [imagr), while data was calibrated 
(using mfeal because the nearby primary calibra- 
tor 1934-638 was used as the secondary calibrator) 
in Miriad. The ATCA does not doppler track, so 



^The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under coopera- 
tive agreement by Associated Universities, Inc 

•^http : //www ■ aips . nrao ■ edu/cook ■ html] 



*The Westerbork Synthesis Radio Telescope is operated by 
the Netherlands Institute for Radio Astronomy ASTRON, 
with support of NWO. 



the correct velocity definition was determined in 
Miriad and applied to the data in AIPS using the 
task cvel. 

We also obtaine d the reduced cubes used 
by iMavbhate et all (|2007l) in their Hi anal- 
reduction strategies are similar, 
Hi emission in the VLA- 
4038, NGC 7252, and 
tails of 



ysis. Their 

optimizing diffuse 

observed tails of NGC 

NGC 3921, and the ATCA-observed 



NGC 3256. Data for NGC 520 w ere obtained 
from iHibbard fc van GorkomI (|1996l ) . A brief de- 



scription of all radio data is presented in Table 1. 
Table 1 includes source distances, beam dimen- 
sions (6i and 62), and beam position angles (P.A.). 
We also list angular and physical pixel areas 
(^pix), channel widths (dv), la detection limits 
in single channel integrated flux densities eg.dv 
(10"* Jy km s~^ channel"^) and column densi- 
ties ejvHi array setups, and original publications 
and/or project codes for the observations. We 
use luminosity distances for all systems (assuming 
Ha = 73 km s^^ kpc"i), except for NGC 4038 
where we adopt the dis t ance measurement of 
Saviane. Hibbard. fc RichI (|2004h . This was done 
for consistency with Mil. All cubes were pre- 
cessed to J2000 for direct comparison with SCC 
positions. 

A cursory glance at Table 1 shows that the pixel 
sizes of the cubes are not a fixed fraction of the 
beam sizes. The size and dimensions of the synthe- 
sized beam depend on the sampling (and weight- 
ing/tapering) of the uv plane in observing these 
sources, which is a function of the array config- 
urations used, observation program design, and 
source declination. Works cited in the Original 
Publications column contain additional informa- 
tion on these issues for each observation. The 
pixel size, meanwhile, is chosen to avoid undersam- 
pling the beam (fewer than «2-3 pixels across the 
beam) and avoid large computational times (^10 
pixels across the beam, depending on the origi- 
nal map size). Different pixel sizes are selected 
between these extremes to maintain sufficient sig- 
nal within pixels and individual channels for the 
varying source intensities of the tidal tail sample. 

3. Analysis 

Our reduced and calibrated cubes contain three 
dimensions: the typical x/y astrometric compo- 



nents, and frequency. Each element of the cube 
registers the intensity (Jy beam~^) at a certain 
position and over a particular channel of observed 
frequency, which tracks the Hi emission over ve- 
locity space. These elements and their errors can 
be transformed into measurements of flux den- 
sity (Jy) by multiplying the intensity by the ratio 
of the pixel and beam angular sizes, Apjx ^bcam' 
whose values are recorded in Tabic 1. 

Our immediate objective is to probe the Hi 
properties on the smallest physical scales afforded 
by our data cubes, specifically where there are 
established cluster candidates. iMavbhate et al.l 



(J2007I ) produced maps of the Hi column density 
for their sample of tidal tails and determined the 
underlying Nfji values for each of their SCC po- 
sitions. Given the differences in resolution scales 
between their optical and radio data, this was in 
essence a pc-to-kpc comparison. Here, we instead 
identify areas in our Hi data that "contain" SCCs 
and compare those to the areas that do not. This 
is a favorable kpc-to-kpc spatial comparison. 

We are therefore concerned with the spectrum 
of Hi emission of different xy positions within the 
smallest regions possible. Individually the resolu- 
tion of each cube is set by the synthesized beam (~ 
several kpc in this sample) , which for the purposes 
of interferometric imaging is sampled by the se- 
lected pixel size (^0.1-2 kpc). In the sections be- 
low we map out Hi characteristics from the cubes 
on a pixel-by-pixel basis. Because this project is 
a synthesis of optical and radio data and termi- 
nology, we will continue to use the term "pixel" 
to describe the xy components of our data cubes. 
Each pixel thus contains a one-dimensional cross 
section of the cube corresponding to all frequency 
(velocity) channels of an individual spatial posi- 
tion. Extracting and characterizing these spectra 
for every pixel with sufficient Hi signal in an es- 
tablished tail region is our critical task, requiring 
careful pixel and channel selection. Not all chan- 
nels of a given pixel will contain Hi signal, which 
necessitates careful masking of "noise" channels 
and selection of "signal channels" to properly as- 
certain Hi properties for that pixel. We describe 
our technique, implemented with the Interactive 
Data Language (IDL), below. 



Table 1. 

Tidal Tail Observations 



— A~~? 

(arcscc ) 



Tail 



Dist.' 

(Mpc 



Ago" 
(Myr) 



Mass 
Ratio^ 



bi' 



Beam 



A ■ ^ 



d-u" 



log 



Project 
Codes' 



Instruments and Arrays 



Original Publications 











The 


MuUan et al. 


(2011) Sample 










NGC 1487E 


10.8 


500 


0.25 


27.0 24.1 16.1 9 

ATCA 375 + 750A + 1.5A + 1.5B + 1.5D + 6A 


0.02 


10 0.93 19.11 
Getts (2001): English & Freeman 


(2001) 


C531 


NGC 1487W 


10.8 


500 


0.25 


27.0 24.1 16.1 9 

ATCA 375 + 750A + 1.5A + 1.5B + 1.5D + 6A 


0.02 


10 0.93 19.11 
Getts (2001): English & Freeman 


(2001) 


C531 


NGC 4747 


20.2 


320 


20 


61.9 48.5 
WSRT 


0.00 


400 


3.84 


16.4 39.47 19.09 
van der Hulst et al. (2001) 




- 


NGC 520 


27.2 


300 


0.05 


24.5 23.8 
VLA C + D 


38.3 


25 


0.43 


10.5 1.13 18.75 
Hibbard & van Gorkom (1996) 


AH412/ 


AH417 


NGC 2992 


36.6 


100 


1 


23.9 17.9 
VLA C 


-14.9 


16 


0.50 


20.5 2.78 19.34 
Due et al. (2000) 




AD402 


NGC 2993 


36.6 


100 


1 


23.9 17.9 
VLA C 


-14.9 


16 


0.50 


20.5 2.78 19.34 
Due et al. (2000) 




AD402 


NGC 2782E 


38.1 


200 


4 


19.8 17.9 
VLA B + C + D 


86.4 


4 


0.14 


10.3 0.53 19.22 AS389/ 
Smith (1997); Smith (1994): Smith (1991) 


AS453 


NGC 2782W 


38.1 


200 


0.25 


19.8 17.9 
VLA B + C + D 


86.4 


4 


0.14 


10.3 0.53 19.22 
Smith (1997): Smith (1994): Smit 


AS389/ 
h (1991) 


AS453 


NGC 2444 


58.2 


100 


0.50 


27.9 27.1 
VLA C + D 


-34.2 


16 


1.28 


20.5 1.70 19.12 
Appleton et al. (1987) 


AA63 


NGC 2535 


59.8 


100 


0.30 


20.4 18.4 
VLA C + D 


16.6 


9 


0.76 


5.1 0.87 19.08 
Kaufman et al. (1997) 




AK327 


NGC 6872E 


62.6 


150 


0.20 


53.5 51.6 25.1 16 
ATCA 750D + EW352 + 1.5G + 750F 


1.46 


19.8 1.64 19.11 
Horellou & Koribalski (2007) 




C979 


NGC 6872W 


62.6 


150 


0.20 


53.5 51.6 25.1 16 
ATCA 750D + EW352 + 1.5G + 750F 


1.46 


19.8 1.64 19.11 
Horellou & Koribalski (2007) 




C979 


Continued on 


next page. . . 





























Table 1. 


Continued 










Tail 


Dist." 
(Mpc) 


Age" 
(Myr) 


Mass 
Ratio"^ 


bi*^ b2'* Beam 
P.A." 


(arc sec"") 


A,„» 
(kpc^) 


dv'^ es^ dv^ 


log 


Project 

Codes^ 




Instruments and Arrays'^ 


Original Publications 











The 


M 


Lilian ct al. (2011) Sample 










NGC 1614N 


65.6 


750 


1 


15.0 10.0 
VLA C + D 




-83.0 9 


0.91 


10.6 3.09 
Hibbard & Yun (1996) 


19.63 


AH0527 


NGC 1614S 


65.6 


750 


1 


15.0 10.0 
VLA C + D 




-83.0 9 


0.91 


10.6 3.09 
Hibbard & Yun (1996) 


19.63 


AH0527 











The Knierman et al. 


(■2003)/Mavbhate et al. 


(2007) Sample 








NGC 4038A 


13.8 


420 




20.7 15.4 
VLA C + D 


24.3 


25 


0.11 


5.2 3.17 
Hibbard ot al. (2001) 


19.20 




AG0516 


NGC 4038B 


13.8 


420 




20.7 15.4 
VLA C + D 


24.3 


25 


0.11 


5.2 3.17 
Hibbard ct al. (2001) 


19.20 




AG0516 


NGC 4038C 


13.8 


420 




20.7 15.4 
VLA C + D 


24.3 


25 


0.11 


5.2 3.17 
Hibbard et al. (2001) 


19.20 




AG0516 


NGC 3256E 


42.8 


400 




25.7 19.3 6.94 
ATCA 1.5B -1- 1.5C -1- 1.5D 


16 


0.69 


33.5 13.06 
Enp|lish et al. (2003) 


20.01 




C061 


NGC 3256W 


42.8 


400 




25.7 19.3 6.94 
ATCA 1.5B -1- 1.5C + 1.5D 


16 


0.69 


33.5 13.06 
Enelish et al. (2003) 


20.01 




C061 


NGC 7252E 


62.2 


730 




26.9 16.1 
VLA C + D 


18.6 


25 


2.28 


42.5 4.53 
Hibbard et al. (1994) 


19.35 


AH0372/ 


AH0412 


NGC 7252W 


62.2 


730 




26.9 16.1 
VLA C 4- D 


18.6 


25 


2.28 


42.5 4.53 
Hibbard et al. (1994) 


19.35 


AH0372/ 


AH0412 


NGC 3921S 


84.5 


460 




19.3 18.1 
VLA C + D 


-85.1 


25 


4.20 


10.7 1.61 18.91 
Hibbard & van Gorkoin (1996) 




AH0417 



Notes. 

"Luminosity distance (Mpc) 

''Dynamical age from Mil and sources therein (Myr) 

'^Defined as M1/M2, with Mi and M2 the mass of the progenitor and perturbing galaxy, respectively 

''Beam dimensions (arcsec) 

"^Beam position angle (deg) 

■'Single pixel area (arcsec^) 

^Single pixel area (kpc^) 

''Channel width (km s~^) 

'la detection limit in the single channel integrated flux density (10^* Jy km s~^ channel"^) 

^Icr detection limit in column density (cm~^ channel"^) 

'^Instruments (ATCA, WSRT, or VLA) and array configurations used in original data acquisition (see original publications for details) 

'ATCA and VLA archive project and proposal codes of acquired data 



3.1. Channel Selection 

We first apply channel selection criteria to all 
pixels, to identify the channels in which Hi emis- 
sion is confidently identified. There are a variety 
of methods to identify signal in individual chan- 
nels of an Hi pixel and to mask channels without. 
For ins tance, THINGS (T he Hi Nearby Galaxy 
Survey; IWalter et al.l 120081 ) required that a signal 
be registered by at least three contiguous chan- 
nels (each of a width <5 km s^^). Their instru- 
ment setup was designed to Nyquist sample the 
warm neutral medium, and they report that this 
method sufhciently highlights real 21-cm emission 
and eliminates noise. 

Our data exhibit a much wider range of chan- 
nel widths, from sa5.13 to 42.5 km s~^, so we must 
avoid such a constraint for cubes where multiple 
channels would span an unrealistically broad ve- 
locity space. Moreover, our data are inhomoge- 
neous in sensitivity, so the appearance of an Hi 
line profile may vary depending on these instru- 
mental constraints and orientation/projection ef- 
fects. Thus, we customized the channel selection 
criteria for each tail after examining the spectra of 
many pixels across the tidal debris. Table 2 lists 
the channel selection requirement as a signal-to- 
noise ratio for single channel signal and equivalent 
integrated flux density for each tail, along with 
the minimum number of adjacent channels A'chan 
needed at that flux density for a pixel to count as 
signal. 

The values in Table 2 were selected after an ex- 
tensive trial-and-error process for each tail that 
produced reasonable spectra and maps of the 
properties described in Section 15^ below. We var- 
ied both Nchan and the necessary signal-to-noise 
ratio until no single, isolated pixels in maps of 
the line-of-sight velocity dispersion (cios) regis- 
tered physically unrealistic values within the tail 
regions (i.e. >150 km s~^, depending on the tail). 
Pixels within a beam size are statistically corre- 
lated, so finding occasional pixels with discrepant 
values in velocity dispersion maps indicated that 
noisy channels outside the velocity space of ac- 
tual emission were unwittingly selected and thus 
tended to increase crios relative to the surrounding 
pixels. The spectra of these locations were also 
visually studied to optimize our masking parame- 
ters. 



Table 2. 

Single Channel Signal Detection Criteria 



Tail 


S/W 


S^uindv' 


iVehan^ 


NGC 1487E 


2.75 


2.56 


2 


NGC 1487W 


2.75 


2.56 


2 


NGC 4747 


3.00 


118.42 


1 


NGC 520 


2.25 


2.55 


3 


NGC 2992 


2.50 


6.95 


2 


NGC 2993 


2.50 


6.95 


2 


NGC 2782E 


2.25 


1.20 


2 


NGC 2782W 


2.25 


1.20 


2 


NGC 2444 


2.00 


3.40 


2 


NGC 2535 


2.25 


1.96 


3 


NGC 6872E 


2.25 


3.70 


2 


NGC 6872W 


2.25 


3.70 


2 


NGC 1614N 


2.25 


6.95 


2 


NGC 1614S 


2.25 


6.95 


2 


NGC 4038A 


2.25 


7.14 


3 


NGC 4038B 


2.75 


8.73 


2 


NGC 4038C 


2.75 


8.73 


2 


NGC 3256E 


3.50 


45.70 


1 


NGC 3256W 


3.50 


45.70 


1 


NGC 7252E 


3.50 


15.86 


1 


NGC 7252W 


3.50 


15.86 


1 


NGC 3921S 


2.00 


3.23 


3 



Notes. 

" Signal- to- noise ratio required for a single channel 
''Minimum integrated flux density required for a single 
channel (10^^ Jy km s^^) 

'^The number of adjacent channels needed at the indicated 
flux to register Hi emission in a pixel 



We then masked out, i.e. ignored, the noise 
channels in the cube, with the exception of 
noise channels alongside channels containing sig- 
nal. These were left unmasked, as were all ad- 
jacent channels until a flux density level consis- 
tent with Jy (± la) was found along either side 
of those original signal channels. This helped to 
prevent clipping of broad spectral features present 
across multiple channels. If any of these additional 
channels registered flux densities below Jy, they 
were then masked to discourage biasing any mea- 
surements toward lower values. When integrated, 
the unmasked pixels provide a total signal typi- 
cally >3ct over the channel-integrated r.m.s. noise. 

We display several examples of results from 
our signal selection procedure in Figure [1] and 



Figure [2] A spectrum from one randomly selected 
pixel from each tail is shown in these figures, plot- 
ted in gray; the overplotted error bars indicate 
the r.m.s. noise. Signal channels are indicated as 
black points. The heterogenous channel widths 
and noisy quality of some of the randomly cho- 
sen Hi profiles are clearly evident in these figures, 
highlighting the need for a varied signal selection 
criteria across the sample. 

3.2. Pixel Selection 

We also limit our analysis to regions in our 
Hi data encompassed by previous optical obser- 
vations, i.e., selecting Hi pixels that are "in-tail" 
and within the WFPC2 FOV of the previous stud- 
ies. Withi n the WFPC2 images o f their datasets, 



Mil and iKnierman et al.l (|2003l ) defined areas 



that were at least one count above the sky level 
in F606W or F555W as "in-tail." These groups 
used the WF2-4 chips in their analysis, choos- 
ing to ignore the planeta ry camera with i t s pro - 



hibitively high readnoise. iMavbhate et al.l ([2007|) 
applied this definition and the optical boundaries 
it demarcates to the Hi maps of their tidal tail 
sample to study the local column densities within. 
In Mil, the expanded tidal tail sample afforded 
a glimpse at tails with different Hi and optical 
morphologies, whose optically bright regions did 
not always coincide with accumulations of Hi (c.f 



Hibbard et al.ll2000l) . Consequently, we redefine 



"in-tail" to refer to an area of the Hi data that 
is confined to the WFPC2 footprint of its corre- 
sponding optical observations, and is either within 
the optical boundaries set by Mil or met the Hi 
criteria of Section 13.11 

For each tail, we determine the physical extent 
of the Hi emission by convolving a map of the in- 
tegrated flux density with an artificial beam twice 
the size of the synthesized beam and masking out 
pixels that do not satisfy the minimum criterion 
for Hi signal. We found that this technique suc- 
cessfully suppresses artifacts in the cube that are 
beyond the true large-scale emission in the tidal 
debris, as opposed to using no smoothing kernel 
or kernels of sizes >3 times that of the beam. Re- 
gions determined to reside at the WFPC2-imaged 
edges of the host galaxies by Mil were rejected. 



3.3. Measured Properties 

3.3.1. Primary Properties 

The three-dimensional nature of the Hi cubes 
provides an opportunity to examine not only the 
distribution of local Hi surface densities, but also 
the velocities and kinematics of the gas. Thus, for 
each in-tail pixel, we compute the total integrated 
intensity /tot (Jy beam~^ km s~^), the mean line- 
of-sight velocity uios (km s~^), and the line-of-sight 
velocity dispersion cios (km s~^). In each pixel, we 
calculate /tot as the zeroth moment of the velocity- 
intensity distribution: 



N 

E 



li dv. 



(1) 



where li is the intensity in each signal channel 
and dv is the channel width in km s^^. Again, 
The corresponding integrated flux density ^tot can 
be obtained by multiplying /tot by the ratio of 
pixel and beam angular sizes Apjx ^bcam (both in 
arcsec^). The line-of-sight velocity and (squared) 
velocity dispersion is given by the first and second 
moments: 



Wlo 






(2) 



and: 



Clo 



1 ^ 



respectively. In Equation (2) and (3), Vi is the 
velocity value for the i-th channel. Summations 
were performed only on channels left unmasked as 
described in Section IXTI13.2I for all equations, and 
errors in li (e/J were propagated through to all 
quantities. We adopt a minimum value of 50% of 
the given channel width for crios (an instrumen- 
tal limit), as well as for errors in both wios and 
CTios • For pixels defined as in-tail by virtue of their 
optical properties alone, these moments were not 
computed and the Hi properties were classified 
as "unfitted." In occasional cases of marginally 
detectable Hi emission (a moment-zero measure- 
ment), higher moments were unable to be reliably 
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Fig. 1.— Randomly selected example spectra of pixels in the Hi debris of NGC 1487E-NGC 6872W. Gray 
points indicate data for all channels; black points indicate channels that have been selected as Hi emission 
with the criteria and procedure of i)3.1l Error bars (r.m.s. noise) are also shown. 



determined, and such pixels were therefore unfit- 
ted in i;ios and uios- 

For visual reference, we display moment-zero 
integrated flux density maps of the full extent of 
each data cube in Figure |3l with WFPC2 fields of 
view overlaid. In all images, north is up, east is to 
the left of the frame, and every WFPC2 footprint 
is w2f47 on a side. The sizes and orientations of 
the synthesized radio beams are represented in the 
bottom left of each frame in this figure. 

We also experimented with finding the equiva- 
lent quantities by fitting the Hi spectra offered by 
each pixel with a series of gaussian functions. In 
practice, we found that our values of Stot for both 
methods agreed very well. However, spectral fit- 
ting often yielded multiple gaussian components, 
so in many cases it could not be easily determined 
which component was "best" associated with the 
tidal debris and which wjos and ctios should be com- 



pared to the moment equivalents (c.f. the com- 
plex line structures evinced in Figures [1] [2]) . The 
intensity- weighted mean of the individual compo- 
nents facilitated a comparison of the wios measure- 
ments, but there was no clear prescription for com- 
paring the number of ctios contributions from the 
gaussian fits to the moment-based measurement. 
Many pixels throughout the tail sample can be well 
fit by multiple gaussian components that overlap 
in varying amounts, so summing the gaussian-fit 
velocity dispersions would be an inaccurate uni- 
versal tracer of kinematic complexity. 

In all, because the results for integrated fiux 
density were very similar, and for simplicity of 
analysis we wish to use single quantities to de- 
scribe the velocity structure at positions along 
a tidal tail, we choose to employ the moment 
measurements for these primary properties. Sim- 
ilar work for galaxy disks can be found in 
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Fig. 2. — Randomly selected example spectra of pixels in the Hi debris of NGC 1614N-NGC 3921S. Gray 
points indicate data for all channels; black points indicate channels that have been selected as Hi emission 
with the criteria and procedure of WS.1\ Error bars (r.m.s. noise) are also shown. 



Tamburro et al.l ( 2009I ). While this method loses 
sight of interesting kinematics in some locations 
along tidal tails, e.g. multiple overlapping gaussian 
fits and extended or asymmetric profiles, complex 
line structures are still reflected in large second 
moments of the intensity distributions. Our mea- 
surement of (Tios is thus still a tracer of such tail 
Hi kinematics. 

3.3.2. Secondary Properties 

Further Hi characteristics can be calculated 
from the moment-based primary properties, e.g., 
the Hi column density, mass surface density, and 
large-scale velocity gradient. The total integrated 
intensity /tot is related to the Hi column den- 
sity Nht (cm~^) by the simpl ified relation (e.g.. 



Verschuur fc Kellermannlll988l ): 



Nhi [cm-2] = 1.247 x 10 



24 



A 



beam 



(4) 



Here, Abeam is the area of the beam in square arc- 
seconds. The error in Njji is given by: 



eA,H,[cm-2] == 1.247 X lO^-* 



/iv: 



Abe 



e/. dv . (5) 



N-a denotes the number of unmasked channels 
for the specified pixel, dv is the channel width 
(km s^^), and e/. is the single channel r.m.s. error 
in li. The corresponding error in Si (esi) is then 
simply {A^i^A^l^J e^. 

We moreover measure the mass surface density, 
kinetic energy density, and velocity gradient of 
each pixel. The mass surface density Smh, and 
its uncertainty csmhi ar e calculated with the foL 



lowing equations (from IVerschuur fc KellermannI 
1988L with some alterations to change masses to 
mass surface densities and include unit conver- 
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NGC 1487E 

«^ /^ NGC 1487W 




NGC 3256W , 



Fig. 3. — Channel-integrated flux density maps of the tidal tails of our sample (see Figure 1 of Mil and 
K03 for optical counterparts). WFPC2 FOVs are overlaid as black footprints; these are all ~ 2f47 on a side. 
The sizes and orientations of the synthesized radio beams are represented in the bottom left of each frame. 
A full resolution version of this figure is available in The Astrophysical Journal. 



sions 



^Mh 



[Mq pc- 



1.00 X lO"* 



At 



and: 



ESmhi [^^O PC" 



1.00x10^ 



^ 



Ay 



(6) 



e/, dv , (7) 



where ej- is the single channel r.m.s. error in in- 
tensity as above. 

To characterize the dynamical state of the Hi 
gas, we further define the kinetic energy density of 
the gas in a single pixel as Ske = ctJ^Mm'^ios'^, 
where the factor a=3 /2 assumes th e velocity dis- 
tribution is isotropic. iTamburro et al. (2009) sim- 
ilarly calculate a pixel-by-pixel kinetic energy for 



their sample of non-interacting galaxies. For 
streaming Hi features at a variety of inclinations 
and projection effects, this constant may differ in 
reality by a factor of a few. Thus we consider this 
only an order-of-magnitude tracer of the actual, 
local turbulent kinetic energy per unit area. 

Lastly, we produce a diagnostic of the locally 
centered line-of-sight velocity gradient across the 
plane of the sky dvios/dr± (km s~^ kpc~^) from 
the root mean square of the velocity gradients be- 
tween the pixel of interest and others in its imme- 
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diate neighborhood. We write this as: 



dvios/dri 



V^ 









p'^ \ i=i 



dr2 



, (8) 



where v is the Hi velocity for the central pixel 
of interest and Vi is the equivalent velocity for 
the i-th nearby pixel with which we construct a 
velocity gradient over dr±^i, the distance in kpc 
between the pixels. A^pix is the number of pix- 
els that surround the central pixel of interest and 
are used in the calculation. Here, we use all pix- 
els with measurable Hi emission that surround 
the central pixel in an annulus whose inner and 
outer radii are set by 0.5 and 1.5 times the maxi- 
mum dimension of the beam. This allows a com- 
parison between the central pixel of interest and 
the closest statistically uncorrelated pixels 1 beam 
away. Thus, while previous quantities pertain to 
pixel scales (^kpc), this velocity gradient is mea- 
sured over several beams, and probes gas motions 
and kinematics on larger (^10 kpc) scales. If wios 
and aios were unable to be determined in a given 
pixel, then that pixel was also deemed "unfitted" 
by dvios/dr±. 

As an example of these calculations, we present 
a montage of Nhi, 5^Mh,) ctios, and wios maps in 
Figure m for the tail NGC 2535. Corresponding 
maps of the remaining tail sample can be found 
in full resolution in the Online Figure Set |4l As 
before, we include the WFPC2 FOV, but also in- 
dicate positions of SCCs with circles. We fur- 
thermore introduce maps of Eke and dvios/dr± 
in this figure. Contours of log Nhi from 20.0 to 
21.4 (cm~^) in steps of 0.2 dex are overlaid in all 
maps. 

4. Relating Hi and SCC Distributions 

Our motivation in this section is to determine 
whether there are observable environmental differ- 
ences between regions of tidal tails that foster star 
cluster formation and those that do not, specifi- 
cally to the best resolution each data set affords. 
For each tail we therefore separated Hi pixels by 
their proclivity for harboring SCCs and compare 
their distributions of properties measured in Sec- 
tion [3l By contrasting the distributions of prop- 
erties of tail pixels that contain no SCCs with tail 
pixels that contain at least one SCC, we assess 



possible environmental differences on a ^kpc scale 
within a given tail. 

However, we realize that not all SCCs are star 
clusters, and some care must be taken in inter- 
preting differences or similarities between H i char- 
acteristics for these two categories. Fortunately, 
if there are quantifiable differences between pix- 
els with SCCs and those without despite contam- 
ination from non-cluster carrying pixels, these are 
still likely substantial indicators of environmental 
dissimilarity. Null results, given this caveat, are 
not as easily understood. 

We reiterate that the pixel size across the sam- 
ple is not a constant fraction of the beam sizes. 
Because we are not yet concerned with rigorously 
comparing the tails (that is deferred to Section[5]), 
this is not an immediate concern. Rather, because 
the resolution remains set by the beam size, vari- 
ations in pixel size for a given tail only affect the 
coarseness of any considered distribution of an Hi 
property and not the distribution itself. Our plots 
of tail properties we discuss later in this section 
and the Appendix use bin sizes that are larger 
than the pixel-to-pixel variations between corre- 
lated pixels, so these are not affected. 

4.1. Individual SCC-Harboring Probabili- 
ties 

In this section, we concentrate on the pixels of 
a given tail that contain at least one SCC. The 
sources detected in Mil are not verified clusters, 
but cluster candidates, so each of these pixels has 
a probability that it contains at least one actual, 
luminous cluster. We therefore first examine the 
source statistics in a manner consistent with Mil 
to determine how likely a source is a cluster rather 
than a background contaminant. From this we 
compute the probability that the pixel containing 
that SCC has at least one cluster. 

Similar to Mil, we determine the relative areas 
within and outside the tail regions surveyed within 
the WFPC2 footprint, ^in and ^out, respectively. 
The main difference in this study, aside from pixel 
sizes, is the optical/Hi tail definition used here. 
Within these areas, we find the number of SCCs 
Nin and Nout- We divide TVjn by Ajn to find the 
overall in-tail source density Ein; these results are 
listed in Table 3 for each tail, along with the in-tail 
completeness limits found in Mil. 
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Fig. 4. — Maps of several derived quantities - Nhi, Smhu ctios, "yios, Ske, and dvios/dr± fit for the tail NGC 
2535. Contours of log Nhi = 20.0-21.4 (cni~^) in steps of 0.2 dex are overlaid on all maps. Although we only 
consider pixels within the WF2-WF4 chips of the overlaid WFPC2 FOV, the rest of the system has been 
mapped for clarity in some cases. Full resolution Figures 4.1-4.22 for all tails are available in the electronic 
edition of The Astrophysical Journal. 
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Mil also calculated the out-of-tail density 
Sout = -^^out/^out and subtracted it from Sin 
to find the background-subtracted, average in- 
tail density of SCCs in these tidal tail regions. 
However, Table 3 shows that many, particularly 
close systems have extensive Hi and optical cov- 
erage and small out-of-tail areas ^out- The num- 
ber of detected out-of-tail objects within these re- 
gions is statistically uncertain, or not at all known 
(^out = 0) in some cases. 

We therefore determine Sout with results from 
star counts genera ted by the Galactic population 



synthesis model of iRobin et all (|2003l l'^l. In a 100 
arcsec^ solid angle centered at coordinates in each 
tail, we generate the number of foreground stars 
visible within the photometric criteria of Section 
[21 This includes aV — I < 1 cutoff and an appar- 
ent y-band magnitude range matching the sources 
detected in, and completeness limits of. Mil (19 
< y < 26). We then determine the resulting num- 
ber of contaminants per out-of-tail area Sout using 
the distance to each tail. We find that variations 
of solid angle (1-100 arcsec^) introduce variations 
of only a few percent. 

In this calculation, we also roughly model the 
photometric errors of the returned catalogs on the 
errors in photometry in Mil. Here, errors in mag- 
nitude m, (Tmag are calculated by a simple lin- 
ear model (Tmag = A + Bm, where A and B are 
0.01 and 0.0035, respectively. Moreover, we em- 
ploy a mean I/-band diffuse absorption of 0.7 mag 
kpc~^, which is recomme nded for intermed iate- 
to-high Galactic latitudes (JRobin et al.l 120031 ). It 
should be noted that this may introduce some un- 
certainty for NGC 3256E/W, which is a partic- 
ularly low-latitude (5 « 11°) object. But ulti- 
mately, it is a more reliable method that is more 
free of statistical uncertainties with low counts of 
detected objects. In addition, these sources were 
screened for extended sources in Mil, so contam- 
ination by background galaxies — already a minor 
issue from the stringent V — I cutoff — is minimal. 
In short, the primary contamination should be 
foreground stars, and this is better assessed with 
the current method. 

We want to concentrate on individual in-tail 
SCCs and their local environments, so we then cal- 
culate the probability Prcai that a SCC-containing 



pixel contains at least one cluster based on the 
probable degree of in-tail contamination by out- 
of-tail sources. To do so we first find the area- 
adjusted number of probable contaminants in the 
in-tail source list, Sout x ^in- This quantity di- 
vided by iVin is simply the fraction of sources 
within the tail that are likely contaminants, or 
Sout/Sin. The corresponding fraction of sources 
that are likely star clusters (SCs) is then /sc = 

1 — (Sout/Sin). 

We record various steps of this process 
in Table 3, along with the calculated F-band 
50% completeness limits for each tail (M 1 1 ; 



Knierman et al.l 120031 : iMavbhate et al.l 120071 ) for 
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reference. For pixels containing one SCC, Preai is 
just the value of /sc for that tail. More generally, 
and especially for pixels with multiple SCCs, Proai 
is the result of applying the binomial distribution, 
using n SCCs in a pixel, each with a probability 
Prcai — /sc- Thus, the odds that a pixel with 
multiple SCCs has at least one actual cluster are 
higher than those listed for single SCC-harboring 
pixels listed in the penultimate column of Tabic 3. 
We indicate in the last column whether Mil found 
strong statistical evidence of tail- wide cluster pop- 
ulation in each of the tail regions. With the ex- 
ception of the low-latitude (and highly foreground 
source-contaminated) NGC 3256W, /sc values are 
all >0.7. In contrast, tails without global cluster 
populations exhibit a wider range of values. 

As a final check, we compared the distribu- 
tion and proper ties of star clusters detected in 
NGC2782W by iTorres-Flores et al.l (|2012[ ) and 
Mil. Because of the instrumental differences be- 
tween the studies (the former use GALEX imag- 
ing and optical spectra, and the latter use HST 
imaging), there is «50% overlap between the 
source catalogs for that tail. That is. Mil de- 
tected about twice as ma ny cluster candidates as 
Torres-Flores et al.l (|2012l ) spectroscopically char- 
acterized clusters. This agrees very well with re- 
sults from the star counts model that predicted 
52% of the Mil detections are contaminants. 

4.2. Comparison of Hi Distributions 

We now compare the distributions of Nhi, 
dvios/dr±, cTjos, and Eke for between Hi pixels 
not containing SCCs and those that do. In doing 
so we must ignore the fraction of pixels that were 
unfitted by the measuring techniques, /unfit; these 
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Table 3. 

Source Statistics 





Source Detections 










Star Counts 






Mil Results 


Tail 


^V,50%" 


N.J 


Naut" 


A,J 


Aout" 


^.J 


Sout» Sout 


xA,J 


/sc' 


Tail-wide 
SCCs?3 


NGC 1487E 


-4.6 


1 





40.60 


5.89 


0.025 


0.039 


1.56 







NGC 1487W 


-4.5 








41.07 


5.52 





0.039 


1.58 


1 




NGC 4747 


-6.0 


6 


1 


84.52 


84.52 


0.071 


0.014 


1.22 


0.80 




NGC 520 


-6.4 


2 





266.21 


28.66 


0.008 


0.007 


1.87 


0.06 




NGC 2992 


-7.2 


10 


3 


305.49 


229.87 


0.033 


0.009 


2.77 


0.72 


/ 


NGC 2993 


-7.3 


11 


1 


367.49 


169.38 


0.030 


0.009 


3.34 


0.70 


/ 


NGC 2782E 


-7.5 


56 


5 


413.93 


165.16 


0.135 


0.006 


2.50 


0.96 


/ 


NGC 2782W 


-7.3 


6 


2 


472.48 


105.25 


0.013 


0.006 


2.85 


0.52 




NGC 2444 


-8.2 


5 


7 


504.38 


850.42 


0.010 


0.004 


1.83 


0.63 




NGC 2535 


-8.4 


36 


8 


871.95 


554.81 


0.041 


0.003 


2.36 


0.93 


/ 


NGC 6872E 


-8.3 


111 


28 


730.59 


824.29 


0.152 


0.009 


6.68 


0.94 


/ 


NGC 6872W 


-8.5 


101 


14 


729.12 


827.22 


0.139 


0.009 


6.66 


0.93 


/ 


NGC 1614N 


-8.6 


19 


8 


304.89 


1408.86 


0.062 


0.003 


1.01 


0.95 


/ 


NGC 1614S 


-8.6 


30 


10 


325.82 


1387.93 


0.092 


0.003 


1.07 


0.96 


/ 


NGC 4038A 


-5.8 


3 





45.79 


30.30 


0.066 


0.050 


2.30 


0.24 




NGC 4038B 


-5.1 


1 





46.91 


29.96 


0.021 


0.050 


2.35 







NGC 4038C 


-5.8 


1 





27.50 


48.71 


0.036 


0.050 


1.38 







NGC 3256E 


-7.8 


32 


9 


460.19 


271.43 


0.070 


0.035 


16.12 


0.50 




NGC 3256W 


-7.8 


51 


4 


580.74 


146.74 


0.088 


0.035 


20.34 


0.60 


/ 


NGC 7252E 


-8.1 


11 


16 


627.03 


928.00 


0.018 


0.003 


2.18 


0.80 




NGC 7252W 


-8.1 


17 


11 


533.54 


1016.92 


0.032 


0.003 


1.85 


0.89 




NGC 3921S 


-8.8 


14 


5 


1718.82 


1147.28 


0.008 


0.001 


1.79 


0.87 





Notes. 

"50% V-band completeness limit (Mil) 
Number of SCCs detected within the tail region 
'^Number of SCC-like objects detected outside the tail region 
■^Area of in-tail region (kpc^) 
'^Area of out-of-tail region (kpc^) 
J In-tail sec surface density (kpc^^) 
90ut-of-tail surface density of SCC-like objects (kpc~^) 
''Probable number of in-tail SCCs that are contaminants 
^Overall fraction of SCCs that are likely star clusters 
■'This column is checked if Mil found evidence for large-scale populations of star clusters within the tail regions 



are typically pixels that fall within the optical tail, 
but have ill-defiiied Hi content. To outline poten- 
tial contrasts between the pixels that are home 
to potential star clusters and those that are not, 
we performed K-S tests between these distribu- 
tions for each measured property. Here, pxs val- 
ues reveal the probability that the SCC-sporting 
and SCC-empty pixels come from the same distri- 
bution of each measurement. If pks ^ 0.01, it is 
likely that these two categories of pixels have two 
separate distributions of an Hi characteristic. The 
results of these tests are also provided in Table 4; 
columns indicate the results of this test performed 
on the indicated property. We indicate where val- 
ues are below 10^^, as specific results below this 
cutoff are not enlightening. 

Overall, it appears that the distributions of our 
calculated Hi column densities, kinematic trac- 



ers, and mechanical energy densities for both 
SCC-carrying and SCC-devoid pixels are indis- 
tinct (pKS > 0.01). NGC 2782E and NGC 6872E 
are remarkable counterexamples, with low com- 
puted K-S probabilities and two of the highest 
tail-wide SCC population densities reported in 
Mil (wO.2 SCCs kpc-2). A few other tails show 
evidence of independent distributions of individ- 
ual properties between these two pixel types (e.g., 
NGC 2992 with Nhi, NGC 2535 and NGC 6872W 
with dvios/dr±, and NGC 3256E and NGC 7252W 
with (Tios).The implications of these results are dis- 
cussed for individual tails in the Appendix. It is 
probable that the varying degree of contamina- 
tion from foreground stars across the tail sample — 
displayed in Table 3 — contributes to a contamina- 
tion of the Hi distributions for the pixels hosting 
at least one SCC and confuses the K-S results. 
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It is additionally important to consider that 
there are minor, inherent shortcomings to K-S 
tests; namely that they are particularly sensitive 
to differences in the centers of cumulative dis- 
tribution functions relative to their edges. Also 
note that many tails have small samples of SCC- 
carrying pixels, to which such statistical tests are 
particularly sensitive, as mentioned at the outset 
of this section, pixel size variations also introduce 
differences in the number of pixels considered, and 
the coarseness of the unbinned Hi measurement 
distributions. This will have some, albeit minor, 
affect on the results of the K-S tests. 

With all uncertainties considered, two distribu- 
tions with K-S probabilities >0.01 may still be 
independent but are difficult to differentiate here. 
Regardless, it appears that the dependence of star 
cluster positions (and possibly clues to their for- 
mation) on individual Hi densities and kinemat- 
ics cannot be easily ascertained on a tail-by-tail, 
property-by-property basis like this. We therefore 
also explore the holistic Hi parameter space of the 
tail pixels, to see if and how all these Hi prop- 
erties collectively relate to SCCs. Figure [5] shows 
the combined distributions of Nhi and crios (left), 
and Eke and dv\os/dr± (right) of NGC 2535. Sim- 
ilar plots for the remaining tail sample are avail- 
able in the Online Figure Set [5] . Bin sizes of 0.1 
dex in log Nhi (cm~^), 5.0 km s^ in uios, 2.0 km 
s~^ kpc~^ in dvios/dr±, and 0.25 dex in log Eke 
(erg pc~^) were employed in contouring. Contour 
levels are scaled to show the fraction of the total 
number of tail pixels examined, including pixels 
that were unfitted by Nhi, cios, dvios/dr±, or Eke 
measurements. We record the fraction of unfitted 
SCC-carrying and SCC-devoid pixels for reference 
as /u^ and /unfit, respectively. Median uncer- 
tainties for all unbinned properties are shown as 
red error bars. 

Dashed lines in the first plot show potentially 
importa nt value s of t he plotted properties — the 
iMavbhate et all (|2007l) Nhi threshold (horizontal 
line), and the turbulent-thermal transition trios ~ 
10 km s~^ (vertical line; this shows where gas ve- 
locity dispersions are driven by turbulent instead 
of thermal motions). The vertical dashed line on 
the second plot indicates a constant dvios/dr± = 
15 km s~^ kpc~^, approximately the value of the 
Oort constant A, the value of the shear in the 
solar neighborhood. This line is purely fiducial 



and meant to guide distinctions between relatively 
shallow and steep velocity gradients. The horizon- 
tal line on that plot shows a Eke cutoff of 10**^ erg 
kpc~^, which considers an energy threshold set by 
setting a hypothetical critical Hi mass surface den- 
sity for star formation of 3 Mq pc~^ (|Sch avc 200J) 
and threshold turbulent velocity dispersion 10 km 
s~^. These demarcations are meant to guide the 
eye and highlight what fraction of a tail has rela- 
tively more dense, turbulent, etc. gas. 

We indicate the fraction of the total number of 
tail pixels that exhibit super- "critical" Hi charac- 
teristics for the combined distributions with /^^ 
for the left plot, and f^ for the right. "Ql" and 
"Q2" refer to particular quadrants delineated by 
the dashed lines, i.e. gas with log Nhi > 20.6 cm~^ 
and (Tios > 10 km s^^ (left plot), and gas with log 
Eke > 46 erg pc^^ and dvios/dr± < 15 km s~^ 
kpc^^ (right plot). Lastly, the SCC-laden pix- 
els are overplotted in both plots as gray points 
if their calculated Proai values are <0.94, and as 
black points if Prcai > 0.94. This value of 0.94 was 
selected with some foresight and will be addressed 
further in Section [5] Recall that Prcai indicates 
the likelihood that an Hi pixel contains at least 
one star cluster. The plots in Online Figure Sets 
4 and 5 are discussed individually for each tail re- 
gion in the Appendix. 

A clear correlation is evident between Hi col- 
umn density and line-of-sight velocity dispersion 
for much of our sample as seen in Online Fig- 
ure Set [S] The example spectra displayed in 
Figures [1] and [2] show that many of these ran- 
domly selected Hi profiles are kinematically com- 
plex, and many would have wide velocity disper- 
sions as measured by the second moment of the 
intensity distribution. The addition of multiple 
Hi components therefore both broadens crios and 
increases Nhi. The degree of overlap between com- 
ponents in a spectrum may depend on the tail 
kinematics and projection effects, so we would per- 
ceive a different relationship — i.e. slope of con- 
toured distributions — between these two proper- 
ties in each tail of Figure Set [5] Most importantly, 
detecting multiple Hi components indicates that 
the actual physical density of Hi at the location of 
a pixel (or SCC) may not be consistently mapped 
by column density in all tails, so some confusion 
in interpreting the importance of Nhi is expected. 
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Fig. 5. — Combined, contoured distribution of measured Hi quantities for NGC 2535, -with levels correspond- 
ing to the fraction of the tail within the indicated parameter space. Pixels with SCCs most likely to have 
star clusters are plotted in black (proai > 0.94), and those less likely in gray (preai < 0.94). The number of 
unfitted pixels, and the tail fraction within the specified quadrants Ql (left) and Q 2 (right) are provided . 
These are defined by the intersection of lines of potential threshold column density ([Mavbhate et al.l 12007 ) 
with the thermal-turbulent velocity dispersion in the left plots; and their combined energy density with the 
value of shear in the solar neighborhood in the right plots. See §4 for details. Full resolution Figures 5.1-5.22 
for all tails are available in the electronic edition of The Astrophysical Journal. 



5. Results and Discussion 

In the previous sections, we analyzed the Hi 
properties to the best possible spatial resolution 
of their individual cubes. We present a thorough 
discussion of each tail's results and figures in the 
Appendix. It is the purpose of that Appendix — 
and the supporting figures presented thus far — to 
review the SCC distributions and Hi character- 
istics of each tail for the benefit of future work 
on individual tidal tails and the interplay between 
their ISM and their star cluster populations. Sev- 
eral intriguing qualitative trends in the ^kpc-scale 
Hi column densities, kinematics, and mechanical 
energy densities can be seen in this sample that 
would not be otherwise ascertained from the sim- 
ple K-S tests of Table 4. These also help direct 
our approach in Section [^3] 

The most obvious trend to note is that sev- 
eral of the systems identified in Mil as hav- 
ing widespread cluster populations (reproduced as 
the last column of Table 3) also tend to have 



a large fraction of their ~kpc-scale Hi with col- 
umn densities exceeding the log Nhi > 20.6 cm^^ 
(NGC 2535, NGC 2782E, NGC 1614N/S, and 
NGC 3256W). This gas is also turbulent {aios 
~ 10-75 km s~^) with velocity dispersions that 
vary tail-to-tail, though probably in part because 
of tail orientation, projection effects, and spatial 
scales. Kinetic energy densities of the tails are 
therefore also generally high (log Eke > 46 erg 
pc~^). The Hi pixels with the highest probabil- 
ity of containing at least one cluster (and usu- 
ally contain more than one SCC) very often lie in 
these high-Nni/EKE regions; the selection of the 
Prcai=0.94 cutoff rcficcts this. Mechanical energy 
density should loosely trace Hi pressures (i.e., in- 
tegrated along the line of sight), so this would im- 
ply that large-scale clustered star formation as de- 
tected by Mil requires high pressures (high Ske) 
and high Nhi on the spatial scales probed here. 
Conversely, a number of tails with low numbers of 
SCCs and low calculated probabilities of their pix- 
els containing real clusters have little high column 
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density and often also little observably turbulent 
gas (NGC 1487E/W if de-projected; NGC 4747, 
NGC 520, and NGC 4038ABC). 

However, many tidal regions with large per- 
centages (^50%) of their Hi exhibiting high col- 
umn densities and kinetic energy densities on these 
~kpc scales do not have pronounced cluster pop- 
ulations in Mil (e.g., NGC 2993, NGC 3256E, 
and NGC 2782W). In some cases tail orientations 
and projection effects may increase the observed 
Nhi, while the actual physical densities of Hi may 
be low. Kinetic energy densities are also high 
from similarly inflated turbulent velocity disper- 
sions (cTios ~ 10-75 km s~^). When compared to 
its eastern counterpart, NGC 2782W presents an 
interesting case for the hidden role of tail- wide H2 
content or dust-to-gas ratios that may also account 
for the observed discrepancy in cluster formation 
efficiencies given similar atomic hydrogen content 
(these two tails can be easily compared because 
they are present in the same Hi cube). The fact 
that NGC 2782E cou ld be the d istorted minor 
galaxy of the merger (Smith 1994^ and h as con- 
firmed CO detections (Smith et al] Il999l ) would 
insinuate that the ISM is very different in these 
two tails. An existing molecular gas reservoir in 
one tail, and/or different dust-to-gas ratios among 
the debris regions would allow effi cient star forma- 



Table 4. 

K— S Probabilities 



tion in one case but not another. Knierman et al 



(|2012n explore this concept further in their assess- 
ment of these tails. 

Moreover, several tails have pixels with high 
probabilities (proai > 0.94) of containing at least 
one cluster in regions of low or undetected Hi 
column density (NGC 4747, NGC 6872E/W, 
NGC 1614N/S). Either Nhi is not important 
for these tails (perhaps in concert with differing 
molecular content), it is important on physical 
scales smaller than we can currently resolve (note 
that these tails have some of the largest angu- 
lar pixel and beam sizes of the sample) , or it has 
been dispersed and/or consumed over ^100 Myr 
timescales, after the clusters had formed. This is 
reasonably consistent with overall Hi column den- 
sity behavior with tail age seen in Mil. Photoion- 
izing feedback from local or centr al sources may 
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also be a factor in these tails (e.g 
I2OOO I. Given our work with cluster probabilities, 
we also cannot say with certainty that many SCCs 
that lie in regions of low Nhi or Ske are not clus- 



Tail 


Nh, 


CTloB 


dvias/<i'>^± 


Ske 


NGC 1487E 




- 




- 


NGC 1487W 


- 


- 


- 


- 


NGC 4747 


- 


- 


- 


- 


NGC 520 


- 






- 


NGC 2992 


1.3C-3 


0.10 


0.52 


0.23 


NGC 2993 


0.17 


0.039 


0.45 


0.086 


NGC 2782E 


<lc-4 


<lc-4 


3.7e-4 


<lo-4 


NGC 2782W 


0.80 


0.43 


0.67 


0.33 


NGC 2444 


0.53 


0.37 


0.090 


0.61 


NGC 2535 


0.012 


0.15 


9.6C-4 


0.039 


NGC 6872E 


<le-4 


<lc-4 


0.076 


<lc-4 


NGC 6872W 


0.24 


0.26 


9.8C-3 


0.33 


NGC 1614N 


- 


- 


- 


- 


NGC 1614S 


0.93 


0.39 


0.017 


0.89 


NGC 4038A 


- 


- 


- 


- 


NGC 4038B 


- 


- 


- 


- 


NGC 4038C 


- 


- 


- 


- 


NGC 3256E 


0.71 


<lc-4 


0.15 


0.61 


NGC 3256W 


0.16 


0.28 


0.011 


0.77 


NGC 7252E 


0.99 


0.50 


0.63 


0.99 


NGC 7252W 


0.42 


7.70-3 


0.044 


0.39 


NGC 3921S 


0.41 


0.27 


0.87 


0.25 



ters (e.g., NGC 7252 and NGC 2444). Feedback, 
interaction dynamics, and gas dispersal may again 
have had an effect here, or these tails may be home 
to smaller-scale star formation events producing 
few bound, luminous clusters as opposed to the 
large-scale triggered events our study is better de- 
signed to address. 

A few tails may also highlight the importance 
of velocity gradients on larger scales. A number 
of tails show large regions of locally, relatively 
shallow dvios/dr± (values vary from the size of 
the beam and the velocity structure of the Hi). 
In several cases these are areas of many differ- 
ent beam physical scales where many SCCs are 
located (NGC 2992, NGC 3256W, NGC 68782E, 
NGC 2535, NGC 2782E, and NGC 7252W), 
or there is a prominent, diffuse optical struc- 
ture that may have lower mass star clusters 
(NGC 1487E/W). This possibly underscores the 
importance of large-scale (~10 kpc) kinematics in 
star cluster formation, including tail-wide shearing 
motions (and lack thereof) that our measurement 
of the line-of-sight velocity gradient might trace in 
certain tails. 

These results are promising, but we again cau- 
tion against more rigorously comparing the Hi dis- 
tributions between the debris regions as shown 
thus far. These tails exhibit a large morpholog- 
ical diversity, both from their range of interaction 
stages and types, and also intractable projection 
effects from inclination and their orientation on 
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the sky. Models show that tails can have a wide 
range of both small- and large-scale thicknesses 
along our line of sight ([Smith 1994); thus sur- 
face densities and other properties discrepant by 
a factor of ^^3 can all correspond to one "in-tail" 
physical density. Highly inclined mergers can also 
dramatically increase the apparent thickness of a 
tail. We see this especially with NGC 1487 and 
NGC 2992, where, for many pixels, their "true" 
Njji values may be lower by an order of magni- 
tude. 

The translation of other line-of-sight proper- 
ties to true physical properties might be even 
more complicated. If the local velocity dispersions 
within tails are not always isotropic (and they are 
certainly not likely to be), the proper interpreta- 
tion of (Tios will vary pixel-to-pixel. Streaming mo- 
tions of multiple Hi components along our line-of- 
sight and beam smearing along the plane of the sky 
will further affect our measurements of crios , not to 
mention dvios/dr±. That quantity bears the addi- 
tional warning that artificial pixel-to-pixel varia- 
tions that arise in trying to define wios in a region 
of low Hi signal (the tail boundaries) will produce 
high dvios/dr± values. This is seen in a few tails in 
Figure Sets |3] and [51 Furthermore, Ske is defined 
as an isotropic energy density, which is a less accu- 
rate description of the actual dynamic state of the 
Hi gas as beam smearing and ordered gas motions 
becomes more prominent. 

Table 3 furthermore reveals that several espe- 
cially distant tails suffer from completeness effects, 
and our source lists for them are incomplete by 
50% or slightly more. In these cases, all pixels have 
a slightly higher probability that they contain an 
sec, with different adjustments for the different 
completeness curves of in-tail and out-of-tail areas. 
The completeness curves of Mil show that out-of- 
tail areas are complete to fainter magnitudes than 
in-tail regions (by ^0.5 mag); thus for a source lu- 
minosity bin close to the completeness limits, we 
are likely to detect more sources out of tail than 
in-tail in distant systems. By this logic our com- 
puted contamination fractions are overestimated, 
and the Prcai probabilities for SCC-carying pixels 
derived from /sc fractions listed in Table 3 are 
too low. Since this effect docs not artificially in- 
crease the probability that pixels with SCCs con- 
tain observable clusters, we consider our existing 
Prcai work to be suitable lower limits. 



Lastly and most importantly, a cursory glance 
at Table 1 reveals that, by necessity, data from 
these tails are quite heterogenous. They arise 
from a number of different observation program 
designs, array configurations and weightings in the 
reduction process, and therefore exhibit a range of 
output channel widths and pixel/beam sizes. Not 
all tails are sensitive to the same velocity disper- 
sions and physical scales, so additional care must 
be taken in assessing gas properties and possible 
thresholds across the tail sample more homoge- 
neously. 

5.1. Comparing Hi, Global SCC Popula- 
tions, and Interaction Characteristics 

In this section we address the issue of hetero- 
geneous physical resolution in our Hi data and fa- 
cilitate a better comparison between the tails in 
our sample. To make our Hi cubes uniform in 
pixel scale, we repeated the work of Section 3- 
4 for the tails, but with their cubes re-gridded to 
the largest physical scale of the tail sample, i.e. the 
2.1 kpc pixel scale of NGC 3921S. We exclude the 
NGC 6872E/W tails in this analysis, as their phys- 
ical beam sizes («16 kpc) are much larger than the 
rest of the sample (several kpc) and probe Hi gas 
on correspondingly different spatial scales. 

For this task, we used the IDL function fre- 
bin, which conserves flux for every spatial slice of 
the cubes. In all calculations requiring the beam 
size, the beam dimensions of each original cube 
was used, with the exception of the calculation 
of dvios/dr± (Equation ([8|)). Here, we used an 
annulus whose inner and outer radii are 4 and 
12 kpc in the determination of iVpix- This simu- 
lates using an 8 kpc beam, which is physically the 
size of the largest beams in the sample (excluding 
NGC 6872E/W) and allowed a more direct com- 
parison across the tail sample. 

We then sorted all the analyzed tail pixels into 
two categories based on pertinent characteristics 
of the interactions — pixels from the list of tails 
with widespread cluster populations as determined 
by Mil vs. those from the tails that did not, pix- 
els from tails that resulted from early-stage and 
late-stage interactions, and pixels from tails pro- 
duced in interactions from different mass ratios of 
their progenitor galaxies. This approach allows 
a more significant assessment of Hi and merger 
properties with respect to cluster populations than 
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additional individual studies would allow, espe- 
cially for the closest half of the tail sample where 
the large physical pixel scales prohibit meaning- 
ful individual study. In the subsections below we 
therefore discuss these distinctions and the distri- 
butions of Hi properties in greater detail, accom- 
panying Figures IMTTl 

In Figures \6\ [3 and [TOl we outline the nor- 
malized distributions of Nhi, dvios/dr±, crios, and 
Eke for pixels from these types of tails. Subcat- 
egories of sec richness, age, and mass ratio are 
represented as blue, filled histograms, and orange, 
outlined histograms. Cumulative distributions of 
these pixels are plotted over the histograms as blue 
and orange curves, respectively. In all cases, bin 
sizes were set by the number of pixels in each tail, 
and the standard deviation of their various distri- 
butions. We also indicate the fraction of pixels 
that were unfitted by the measuring techniques as 
/unfit; as before these are typically pixels that fall 
within the optical tail, but have ill-defined Hi con- 
tent. 

We also perform K-S tests between the (un- 
binned) blue and orange distributions in each 
panel of these plots. In all cases we calculate 
extremely low probabilities {pks < 10~^) that 
these pairs of distributions are statistically indis- 
tinguishable. Thus, the differences in Hi proper- 
ties between tail types we note below are likely 
real and significant. 

Figures [71 m and [11] are meanwhile analogous to 
the contour plots of Online Figure Set [5] excluding 
the presentation of fsc, which changes tail-to-tail. 
Individual pixel Preai probabilities were computed 
as before for each pixel containing at least one 
sec from each tail and are reflected here. The 
fraction of unfitted SCC-carrying pixels /y^^ is 
also shown. Contour bins are also the same as in 
Figure [5l as are all definitions of Hi thresholds. In 
each of these figures, we juxtapose the Nhi-ctios 
and 'EKE-dvios / dr ± parameter space of the two 
types of tails considered for visual comparison. 

In the following subsections, it is important to 
note that projection effects are still an issue. How- 
ever, a reasonably random set of tail orientations 
would artificially change values of properties by 
different amounts and blur distinctions between 
distributions of an Hi variable between two types 
of tails. We therefore contend that, while statis- 
tically similar distributions may be inconclusive 



about physical differences between types of debris, 
statistically different distributions are still signif- 
icant even against irrevocable geometric effects. 
Any differences between tail types are likely telling 
of real physical dissimilarities in the Hi content 
between them. 

5.1.1. Hi and Tail-wide SCCs 



Of the Mll/ lKnierman et all (|2003l ) tail sam- 
ple with Hi coverage. Mil identified signifi- 
cant tail-wide cluster populations (to 2.5cr con- 
fidence) in NGC 2993, NGC 2782E, NGC 2535, 
NGC 6872E/W, NGC1614N/S, and NGC 3256W. 
While NGC 2992 did not have a strictly global 
sec population, many of its cluster candidates are 
found towards the tip of the optical tail, implying 
(along with their photometry and present statis- 
tics) that many could be real clusters. Therefore, 
we include that system along with the rest of those 
listed as tails with very likely widespread cluster 
populations. We therefore compared the Hi prop- 
erties of the re-gridded pixels of these tails with 
the remaining, globally SCC-poor, tail sample as 
in Section l4?2l 

We present the results of this investigation in 
Figure [6l with solid blue and empty orange nor- 
malized histograms representing pixels from SCC- 
poor and SCC-rich tails, respectively. Most strik- 
ingly, the distributions of Nhi and Eke are very 
clearly different between putatively cluster-rich 
and cluster-poor tails — the former peaks at about 
0.5 dex higher in Nhi than the latter, with large 
percentage of its measurable Hi pixels above the 
Maybhate et al.l {2007) threshold. Tails with large 
numbers of SCCs also on average have higher Hi 
kinetic energy densities on ~kpc scales (again, by 
~0.5 dex). Despite individual cases of localized, 
probable clusters in tails with low column den- 
sities (ss20% of the pixels of SCC-rich tails are 
unfitted for individual Hi properties, compared 
to «10% for SCC-poor tails), Nhi and Eke are 
still apparently important variables in the holis- 
tic appraisal of tidal tail cluster formation ac- 
tivity. These variables could be tied to physi- 
cal Hi densities and pressures, which likely play 
prominent roles in the formatio n of star clusters 
and other bound structures ( e.g., Elmegreenlll993t 
Elmegreen k Efremovlll997tlElmegreenll2008[ ). 



Similarly, SCC-rich tails have an extended dis- 
tribution of CTios , with many pixels showing values 
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Fig. 6. — Normalized histograms of Nhi, dvios/dr±, crios, and Ske for identically-sized pixels from tails that 
do not have global SCCs populations (blue) and those that do (orange). The fraction of unfitted pixels /unfit 
is also recorded (see ^S.l.ll for details). A full resolution version of this figure is available in The Astrophysical 
Journal. 
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in the 30-75 km s~^ range. These tails seem to 
have more highly disturbed gas. While relatively 
more Hi in SCC-poor tails appears to have line- 
of-sight velocity gradients <10 km s~^ kpc~^, the 
velocity gradient measurement distribution trails 
off to generally higher values than that of the SCC- 
rich tails. This is also evident in contrasting the 
contour plots of Figure [T] the top row displays 
the contoured distributions for SCC-rich tails, and 
the bottom row the equivalent plots for SCC-poor 
tails. While about 60% of SCC-rich tail Hi is mea- 
surably turbulent and of high column density, the 
percentage for the SCC-poor tails is only about 
20%. 

Understandably, SCC-rich tails have more over- 
plotted points than SCC-poor tails, and especially 
more pixels containing high probability clusters. 
Again, many of the highest probability clusters lie 
in regions of high column and kinetic energy den- 
sity. However, as noted in the Appendix, many 
tails have probable clusters in unfitted or Iow-Nhi 
pixels. Given the statistics of many of the sources, 
it is likely that many of the gray points also indi- 
cate positions of real clusters in this Hi parameter 
space. Evidently, these plots provide a compelling 
argument for the importance of kpc-scale Hi in 
tidal tail star cluster formation as a whole, but 
there are individual exceptions to the rule. 

5.1.2. Hi and Tail Dynamical Age 

There is much evidence that tail age or inter- 
action stage has a pronounced effect on star for- 
mation phenomenolog y. In their sara ple of com- 
pact groups, Martinc z-Badenes et al.l ([2012) find 
evidence for a general decrease of molecular gas 
deficiency with the evolutionary phase (age) of 
the interacting group. They interpret this as an 
enhanced conversion of atomic to molecular gas 
early in the interactions, followed by gradual strip- 
ping and dispersal of Hi in their later phases. 
In their models of compressive tide s in pairwise 
galaxy mergers, iRenaud et al.l (|2009l ) additionally 
show that, while regions with tidal compression 
can appear in interactions at multiple stages, they 
often occur at greater frequency in younger stages. 
Other models of star-forming tidal debris often in- 
dicate that strong bursts of star formation occur 
early in interactions, soon a f ter periap se (IChien 



20101: iDi Matteo et al.l l2008l iMihos fc Hcrnauist 



19941 ). Perhaps as a consequence. Mil noted that 



the most successfully star cluster-forming tails of 
their sample were all relatively young (^250 Myr), 
although they did not find a statistically signifi- 
cant difference in the age distributions of SCC-rich 
and SCC-poor tails. 

In Figure [8] we exhibit normalized histograms 
of Hi properties of identically sized pixels from 
young (<250 Myr; blue) and old (>250 Myr; or- 
ange) tails. Ages are from Mil and sources therein 
(see Table 1); the age cut was selected to ap- 
proximate the median of the Mil age distribution 
(«300 Myr). For all examined properties, these 
two categories of tails are typically quite differ- 
ent; younger tails enjoy higher column densities, 
have less gas with high velocity gradients (>15 
km s~^ kpc~^) over relatively large scales (~10 
kpc), but have a wider range of local (^kpc-scale) 
velocity dispersions. There is more gas with a\os 
measurements >20 km s~^. The fractions of unfit- 
ted pixels for these Hi characteristics tend also to 
be lower in younger tails (/unfit ~ 0.06, compared 
to «0. 15-0. 23 for older tails). Consequently, the 
Hi of younger tails have higher mechanical energy 
densities on average. In all, the Hi gas in younger 
tails may be under higher pressures and is more 
kinematically disturbed than that of older tails. 

Figure |9] contrasts the contoured Hi parameter 
space of young and old tails in the top and bottom 
rows, respectively, along with their SCC-carying 
pixels. There are more probable SCCs in younger 
tails, despite the larger number of older tails in 
our sample (by a factor of sa2). Their host pixels 
follow the higher percentage of higher column den- 
sity and turbulent gas (sa52% compared to 27%), 
as well as the higher percentage of high mechan- 
ical energy densities and low large-scale velocity 
gradients («86% compared to 43%). As before, 
there are pixels with possible SCCs in Iow-Nhi ar- 
eas in tails of all ages; the fraction of these with 
unfitted Hi characteristics is higher in older tails 
than in younger tails. 

Debris like NGC 1614N/S are examples of older 
tidal regions with ongoing or prolonged star clus- 
ter formation, broadening the distributions of Hi 
properties of older tails and contributing high- 
probability cluster-carrying pixels in the bottom 
plots of Figure m The unplotted tail NGC 6872W 
shows most of the black points in the "subcritical" 
regions of the left plot of its contour diagram in the 
Online Figure Set [5l perhaps as a consequence of 
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its large beam size (diluting smaller-scale Nhi con- 
centrations) or rapid Hi depletion. Both of these 
cases are described in the Appendix. But in gen- 
eral these results are consistent with findings in the 
literature reporting the significance of young ages 
in establishing the necessary conditions — ambient 
densities and pressures — for star forming behav- 
ior. 

Shearing kinematics in tidal debris also help ex- 
plain the dispersion of SCC positions in especially 
older tails. Given the witnessed velocity gradi- 
ents of galactic disks and debris regions examined 
here («0-20 km s^^ kpc^^ on average), it is likely 
that star clusters migrate a significant distance 
from their "original" Hi pixels by the time they 
are observed. An order-of-magnitude estimate re- 
veals that clusters can move several pixels (several 
kpc) away from their original pixel by the cutoff 
age of 250 Myr for these plots. Combined with 
age-dependent Hi depletion and the tendency of 
clusters to fade as their stellar populations evolve 
(e.g., Mil and references therein), a tighter cor- 
relation between SCC detections and Hi column 
densities and kinematics should logically be ex- 
pected for young tails than older tails. 

5.1.3. Hi and Progenitor Mass Ratio 

Major mergers outnumber minor mergers in our 
sample by a factor of ~2. The normalized his- 
tograms of Figure [TU] show that tails from minor 
mergers seem to have a larger fraction of relatively 
high column density Hi than do tails from major 
mergers. Most of their gas, however, has line-of- 
sight velocity dispersions <40 km s~^, while major 
mergers have a slightly more extended distribution 
of (Tios out to sa50 km s^^. The overall appearance 
of the CTios histograms for these two types of tails 
is somewhat similar, but appears statistically in- 
dependent. Given the higher fraction of Hi with 
column densities >10^°'^ cm~^ in our subsample 
of tails from minor encounters, this tail type also 
enjoys a bump in kinetic energy density histogram 
towards log Ske ^ 46.5-47 erg pc~^. 

The distributions of dv\os/dr±_ are the most 
strikingly different between tails from major 
and minor encounters. Velocity gradients ap- 
pear largely confined to values below «10 
km s~^ kpc~^ on wlO kpc scales for the latter cat- 
egory, while our set of tails from major mergers 
has a drastically broader range oidv\os/ dr±^ values. 



Furthermore, tails from major mergers have a con- 
sistently higher fraction of unfitted pixels («20%) 
than tails from minor mergers («2%). 

Many of these trends are of course refiected in 
the associated contoured distributions of these two 
types of tails in Figure [TTl with the addition of 
SCC-carrying pixels as before. There are clearly 
more cluster candidates of both relatively high and 
low probability in the tails of major interactions 
than in minor ones, but again tails from major 
interactions outnumber those from minor interac- 
tions 14 to 6 in this figure. Despite much of the Hi 
gas inhabiting the ostensibly subcritical column 
density region of the parameter space explored in 
the top row (major encounters) of Figure llll many 
of the pixels with high probabilities of having at 
least one cluster still lie in the higher column den- 
sity region. 

Between the evident trend of minor merger tails 
having higher Njji overall and/or a smaller range 
of dvios/dr± for this medium, it is possible that 
major mergers are prone to more intense interac- 
tions. A higher Hi dispersal efficiency may be pos- 
sible in major interactions; the larger fractions of 
unfitted pixels in tails from these encounters might 
testify to this. Stronger interactions may also en- 
gender greater disruptions to the velocity fields 
on ~10 kpc scales. For example, iKeel fc Bornd 



(|2003[ ) find that the major merger NGC 6621/2 
exhibits strong Ha velocity perturbations and has 
an extensive population of star clusters. This con- 
trasts with their observations of NGC 5752/4, a 
minor merger with evidently less kinematic dis- 
ruption. While the interaction ages of these two 
collisions are different (~100 Myr for the fo rmer 
and ^250 Myr for the latter), iKeel fc Bornd con- 
tend that star cluster formation depends on a mass 
ratio-modulated strength of the interaction more 
than the interaction stage. However, they also 
examine cluster populations across the progenitor 
galaxy disks rather than just the tidal debris, mak- 
ing these studies difficult to compare. 

But most critically, it must be noted that the 
single and combined distributions of Hi properties 
of tail pixels from minor and major mergers bear 
a qualitative resemblance to those from young and 
old tails, respectively. This can be seen by compar- 
ing Figure |9] with Figure [TTl Table 1 shows that 
this results from the limited sample size; while 
the Mil sample is reasonably mixed in terms of 
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Fig. 10. — Normalized histograms of Nhi, dv\os/dr±, cios, and Ske for identically-sized pixels from tails 
from minor (M2/M1 < 1; blue) and major {M2/ Mi > 1; orange) encounters in our sample. The fraction of 
unfitted pixels /unfit is also recorded (see ^5.1.3l for details). A full resolution version of this figure is available 
in The Astrophysical Journal. 
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age and mass ratio, th e subsample attributed to 



Knierman et alJ ([20031) supplies a large number of 



older tails from major encounters. Thus, while our 
observations and analyses can be viewed through 
the lens of interaction age and/or mass ratio, our 
finite sample size makes it difficult to ascertain 
the effects of one or the other. Removing the 
iKnierman et al.l ( 20031 ) sample for this work would 
only further decrease the sample size and make 
any results more sensitive to particular tails and 
their unstudied additional effects that certainly 
contribute to their SCC populations. We there- 
fore opt not to do that here. Rather, we maintain 
that a combination of tail age and progenitor mass 
ratio (and likely a number of other variables not 
studied here) is at least partly responsible for the 
Hi and SCC phenomenology we present. 

5.2. Turbulence: A Cause or Consequence 
of Star Formation? 

Supernovae (SNe), particularly Type II, are 
generally regarded as the dominant source 
of fe edback in the ISM in galactic contexts 
fe.g.. Ivan den Bosc hll2000l: iTasker fc BrvanlbOOSi : 
Tasker fc Brvanl 12006: Spitzer 199Q|). It then be- 



comes a question for this project of whether the su- 
pernovae produced in prior star formation events 
are the primary source of the mechanical ener- 
gies observed in the neutral Hi in these tidal tails. 
That is, do the kinetic energy densities (and pre- 
sumably pressures) we measure result from star 
and cluster formation, or are they the cause of 
star and cluster formation? 

Using UV/IR-derived star formation rates 
(SFRs) of a sample of relatively low incli- 
nation disk ga l axies in the THINGS survey, 
iTamburro et al.l (J2009I ) calculated the pixel-by- 
pixel kinetic energy densities that would be pro- 
duced from feedback from core-collapse SNe. In 
comparing these to measurements of the actual 
Hi energy densities, they find that the efficiency 
of energy injection from SNe would have to be un- 
realistically high in order for their measured Ske 
values to be supplied by star formation. Not hav- 
ing a uniform survey of UV, IR, or Ha diagnos- 
tics of star formation rates for our sample, we in- 
stead aim to determine the potential star forma- 
tion rates of our tidal tails from their Eke and 
suggest whether the results are reasonable for this 
kind of environment. 



As in iTamburro et al.l ( 2009f) . we begin by 
equating the kinetic energy density to that the- 
oretically provided by supernovae; i.e. Ske — 
esN-E'SNT'SN- ^-SN is the rate at which energy is 
released by SNe, which is a function of the SFR; 
esN is the efficiency with which that energy is 
converted to the observed turbulence, and tsn is 
the turbulence dissipation rate. Combining their 
Equations (5)--(7), characterizing the dissipation 
timescale and fraction of the stellar populations 
that produce SNe, we derive the following expres- 
sion for the pixel- by-pixel SFR (T): 



T 



Mq yr-i 



2.03 X 10" 



eSN(l + /la) 



(a^ 



SN 



\{M/Mq[ 

Mhi 



A 



100 pc 



Ma 



O^lo 



10 km s 



(9) 



As before, a corresponds to the degree of 
isotropy (here assumed to be 3/2), and the fac- 
tor l/(l-f /la) is a term included in our derivation 
to account for the fractional contribution of type- 
la SNe. Following iMannucci et al.l (|2005l ). we con- 
sider /la values in the range 0-0.33. The term /sn 
indicates the fraction of (recently formed) stars 
that end as core-collapse SNe, while (M/Mq) de- 
notes the average mass of the stellar population 
within the examined region. Dividing these two 
terms in the equation, we consider values in the 
range 0.9 x 10 ~ ^ to 1.3 x 10^^, as shown by 
Tamburro et al.l ( 20091 ) for an upper IMF mass 
limit of 20-50 Mq, with a mass function of the 
form (j){m) = to"""; ao = 1.3 for 0.1-0.5 Mp, and 
2.3 for 0.5-120 M© (e.g.. ICalzetti et al.ll2007l and 
references therein). We use values in the range 
0.1-0.5 for the efficiency esN] the l ower limit is sug- 



geste d by numerical simulations (jThornton et al 



19981 ). and the upper limit provides a conser- 
vative lower limit to T. The turbulence driv- 
ing scale is given by A , which simulations pre- 
dict to be 100±30 PC (jjoung fc Mac Lowl [200I 



de Avillez fc Breitschwerdtll2007l ). The mass of Hi 
is Mhi, and trios as defined here. Note the addi- 
tional factor of a\os', this results from the implicit 
expression of tsn in terms of A and ctios. The con- 
stant 2.03 X 10""^ incorporates all unit conversions 
and the 10^^ erg that are as sumed to be released 
from each type II SN event (JHeileslll987[ ). 
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Table 5. 

othctical Star Formation Rates 





Variable P 


IX. Scale 


Const. P 


X. Scale 


Tail 


mill. T 


max. T 


min. T 


max. T 




(Afo 


(Mq 


(Mo 


(Mq 




yr-') 


yr-^) 


yr-i) 


yr-i) 


NGC 1487E 


0.056 


1.0 


0.15 


2.7 


NGC 1487W 


0.035 


0.62 


0.069 


1.2 


NGC 4747 


0.0079 


0.14 


0.0053 


0.095 


NGC 520 


0.080 


1.4 


0.080 


1.4 


NGC 2992 


9.1 


160 


11. 


190 


NGC 2993 


3.0 


53. 


2.8 


50. 


NGC 2782E 


2.6 


47. 


3.1 


56. 


NGC 2782W 


0.53 


9.4 


0.54 


9.7 


NGC 2444 


1.1 


20. 


1.3 


23. 


NGC 2535 


1.5 


28. 


1.9 


33. 


NGC 6872E 


18. 


320 


18. 


320 


NGC 6872W 


15. 


260 


14. 


260 


NGC 1614N 


0.44 


7.9 


0.48 


8.6 


NGC 1614S 


0.24 


4.3 


0.48 


8.6 


NGC 4038A 


0.022 


0.39 


0.0082 


0.15 


NGC 4038B 


0.0011 


0.020 


0.0019 


0.033 


NGC 4038C 


2.2e-4 


0.0039 


2.4e-4 


0.0042 


NGC 3256E 


2.3 


41. 


3.2 


57. 


NGC 3256W 


6.0 


110 


7.4 


130 


NGC 7252E 


0.26 


4.6 


0.26 


4.6 


NGC 7252W 


2.4 


43. 


2.9 


52. 


NGC 3921S 


1.5 


26. 


1.5 


26. 



We then mapped out the pixel-by-pixel SFR 
of all the tails using our maps of Smh, and cios, 
for the combinations of aforementioned constants 
that yield the smallest and greatest result. We 
performed this for both original maps and those 
produced at the constant pixel scale. In prac- 
tice, changing the efficiency of injection esN has 
the dominant affect in this equation; we added 
the results of all pixels belonging to individual 
tails and record the results in Table 5. We posit 
that if the tail-wide ongoing star formation rate 
(of similar timescale as the turbulence dissipation 
timescale) could fall between the ranges quoted in 
this table, it is possible that star formation and the 
SNe explosions that result from it generally caused 
the turbulent kinetic energy densities we observe 
rather than was an effect of it. Ideally, uniformly 
deep, wide-field imaging in Ha of our tail sample 
would provide the star formation rate on the corre- 
sponding timescale (^10 Myr), allowing a direct 
and straightforward comparison. Unfortunately, 
no such survey exists. 

There are a fe w potentially elucidati ng case 



studies, however. iKnierman et alj (|2012l ) find a 



Ha-derived SFR for the entire western tail of 
NGC 2782 of 9 X IQ-^ Mq yr"! kpc'^. Adjust- 
ing for the in-tail area used here, we calculate a 
rate T «2 x 10"'^ Mq yr~^. This is several or- 
ders of magnitude lower than the values reported 



in Table 5; thus the Hi turbulence is not likely 
caused by star form ation in this case. Moreover, 
the Ha data of .Mihos et alJ ([1993) indicate that 
most current star formation in NGC 6872 is con- 
fined to its tidal tails, and is occurring at a rate 
~3 Mq yr~^. The individual rates for the eastern 
and western tails must be smaller; but in either 
case the rate is far below the minimum T needed 
(~15 Mq yr~^ for either tail and pixel size) for 
the observed turbulent energy to be supplied by 
feedback from star formation overall. 

While the turbulence we observe is purely on 
~kpc scales, it is likely crucial to the kinematic 
nature of the ISM on smaller, star-forming (^pc) 
scales. Turbulence is injected on hundreds of pc 
to ~ k pc-scales (e.g., Mac Low! Il999t Elmegreen 

cascades down to 



l1 



1993; iBlock et al 



2ttlCi) . and 
smaller scales until it is dissipated into ther 



mal energy by vi scosity ( Elmegreen et al.l l2003t 



Kolmogorovl 119411) . There are a number of ob- 



served radiative pathways available to eliminate 
this energy effectively and faster than the turbu- 
lence decay timescale by several orders of mag- 
nitude; e.g., collis ionally-excited lines of Civ 



( Wolfire et al.l 20031). rovi brational transitions of 



H2 (eg. iCluver et al. 201ff), and radio/X-ray con- 
tinuum emission (Kaufman et al. 2012). 

In effect, we are left with a top-down evolu- 
tionary process, where gas kinematics at observ- 
ably large scales — likely influenced by a combi- 
nation of tail age and progenitor mass ratio as 
seen in Section 15.11 — influences star formation on 
the otherwise unresolved, smaller scales. Tur- 
bulent gas generates a log-normal density prob- 
ability densi t y fun ction (PDF) in the local ISM 



(|Wada et al.l I2002I ). whose width — dispersion in 
increases mainly with turbul ent Mach 

Thus, 



density 

number ( Krumholz &: ThompsonI 120071 ). 
turbulence on large scales affects the density distri- 
bution of gas on small scales. This is exacerbated 
in mergers, where additional, expeditious kine- 
matic perturbations create deviations in the log- 
nor mal distribution, favoring more high-density 
gas (JBournaudll201lh . 

In all, our results present a compelling case 
for the importance of interaction-established, kpc- 
scale Hi kinematics and column densities in deter- 
mining the capacity of pc-scale star cluster forma- 
tion in tidal debris. Other, currently unsurveyed 
interaction properties — e.g., dark matter halo po- 
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tentials, interaction speeds and impact parame- 
ters, H2 content, and dust-to-gas ratios — may help 
set the stage for luminous star cluster formation 
in tidal tails and may help account for the vari- 
ety and distribution of SCC populations witnessed 
here. However, the character of the neutral hydro- 
gen medium, described by its densities, pressures, 
and kinematics, also appears to be an invaluable 
factor in this process. We have shown that high- 
confidence SCC detection often occurs in tails with 
large fractions of high column density, turbulent 
gas (Section 15 .1.11) . Both the putative star cluster 
populations and the underlying Hi medium clearly 
varies depending on the age (Section l5.1.2|) and/or 
the progenitor galaxy mass ratio (Section 15.1.31 ) 
of the interaction, indicating that these dynamical 
properties are important in tidal tail star cluster 
formation as well. 

6. Conclusions 

In this paper, we have performed a pixel-by- 
pixel characterization and analysis of the Hi con- 
tent of 22 tidal debris fields from pairwise galaxy 
interactions, comparing these properties between 
pixels that contain star cluster candidates (SCCs) 
and those that do not. In our tail sample we mea- 
sured Hi column density Nhi, mass density ^Mh^ 
line-of-sight velocities wios and velocity dispersions 
(Tios, as well as mechanical energy densities Eke to 
the best pixel-scale spatial resolutions (^kpc) our 
archival data afford. We also constructed a tracer 
of the locally-centered line-of-sight velocity gradi- 
ent dvios/dr± for these regions, examining the ve- 
locity fields of the tails over larger, ^10 kpc scales. 

While tail gas kinematics and geometry, small- 
scale star and cluster formation events, and feed- 
back likely produce dispersion in our results, it 
appears that many "successful" kpc-scale trig- 
gered clustered star formation events often take 
place in ~kpc-scale regions of high Hi gas col- 
umn and kinetic energy densities. Kinetic energy 
density should generally trace Hi pressures (inte- 
grated along the line of sight), so this would im- 
ply that clustered star formation as detected by 
Mil requires high pressures (high Ske from mod- 
erately turbulent gas), but also high Nhi on these 
scales. It is possible that these conditions could 
arise from intermittent, compressive tides (e.g.. 



Renaud et al. 2009 ) or gravitational instabilities 
(;e.g.. lBournaudll2010l ). 

We repeated our analysis for data regridded to 
a constant pixel scale, and compared the distribu- 
tions of Hi properties between tails of different 
global cluster populations, dynamical ages, and 
progenitor mass ratios. Despite variation in indi- 
vidual tail phenomenology and SCC populations, 
several clear trends and contrasts are evident in 
the Hi content of these types of tails: 

1. Tails with global cluster populations (as 
identified in Mil) tend to have Hi gas of 
higher column and kinetic energy densities 
than tails with no apparent tail-wide cluster 
populations, with much of it (>50%) above 
the putative log Nhi = 20.6 cm~^, log Eke = 
46 erg pc~^ thresholds mentioned here. The 
distribution of velocity dispersions also has 
a greater contribution of high values (trios ^ 
30 km s^i). 

2. Relatively young tails (dynamical age <250 
Myr) have more pixels with high probabil- 
ity clusters than older tails (>250 Myr), and 
contain gas with larger column and kinetic 
energy densities on average. There is also 
more gas overall at high velocity dispersions 
(c^ios ^ 30 km s~^) in these types of tails. 
Older tails are more likely to be affected by 
age-dependent fading in their cluster popu- 
lations as well. 

3. Tails from minor encounters appear to have 
more relatively high column density Hi than 
tails from major interactions. When mea- 
sured over a common, ^10 kpc-scale veloc- 
ity field, major encounters appear typically 
slightly more capable of driving the gas to 
higher (>10 km s~^ kpc^^) velocity gradi- 
ents. The strength of the encounter, gauged 
in part by the interacting mass ratio (among 
other variables) may help control the effi- 
ciency of Hi dispersal and disruption of the 
velocity fields in its resulting tidal debris. 
Major mergers may drive gas to high kpc- 
scale velocity dispersions, but may also en- 
courage stronger local depletion and larger- 
scale shearing motions in certain cases. 

4. Because of the strong degeneracy in our sam- 
ple between old tails and tails from major 
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interactions, it cannot be well determined 
whether the observations noted above are 
explained by primarily age or mass ratio, or 
a combination of the two. Other unstudied 
properties like interaction speed, dark mat- 
ter halo structure, and ISM dust-to-gas ra- 
tios are likely to also play roles in star cluster 
formation in tidal tails. 

5. In calculating the required star formation 
rate (SFR) for the turbulent energies to be 
provided by feedback from star formation, 
we deduce that it is unlikely that the kinetic 
energy densities of the Hi medium are a con- 
sequence of star formation, but rather a pre- 
existing environmental condition. However, 
we can only verify this for a few systems that 
have data of SFR indicators (i.e. Ha) that 
trace activity around the same timeframe as 
the ~10 Myr turbulence injection and dissi- 
pation. 

While these results are promising, much work 
remains in resolving certain issues. Additional op- 
tical photometry (e.g., broadband UB coverage to 
add to existing VI images) is required to accu- 
rately age-date cluster candidates for a more pre- 
cise comparison to Hi properties than the prob- 
abilistic work incorporated here. A larger tail 
sample would also allow a more robust compar- 
ison between Hi measurements between certain 
tail types, and additional subcategories within. A 
homogeneous Ha survey of the tidal debris sam- 
ple would facilitate a quantitative comparison be- 
tween hypothetical star formation rates calculated 
here, and their measured equivalents over the ap- 
propriate «10 Myr timescale. Lastly, consistent 
modeling of all interacting systems would provide 
a direct link from the projected quantities we ob- 
serve here to real physical properties like Hi num- 
ber density and local pressure. 
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A. Notes on Individual Tails 

Below we present an overview of the observed 
Hi properties of eaclr tail as evinced in Figure |3] 
and online Figure Sets 2] and [5] For more infor- 
mation about these systems and their SCCs, con- 
sult Mil and references therein; several dynamical 
and observed properties from those works are oc- 
casionally mentioned here. 

A.l. NGC 1487E 

This eastern tail of a «500 Myr old 
merger remnant b etween similarly massed galaxies 
(|Lee fc Lee|[2005l ). NGC 1487E contains a single 
sec that is not likely real (/sc = 0). This system 
is the closest of the tail sample, so optical images 
from Mil reveal that star formation may have oc- 
curred in detectable, fainter structures. Some of 
these may be ~10^ Mq star clusters formed in 
situ. This tail clearly favors the formation of dis- 
persed, lower-mass structures over more massive 
star clusters. 

Figures 4.1 and 5.1 display a high column den- 
sity Hi tail on ~kpc scales, with line-of-sight ve- 
locity dispersions typically ^20 km s^^. The ve- 
locity field of this tail varies smoothly and pro- 
duces line-of-sight velocity gradients within this 
tail all >5 km s^^ kpc^^. Figure 5.1 indicates that 
NGC 1487E has a high fraction of turbulent and 
dense gas, but Mil suspected that inclination ef- 
fects artificially increase the apparent Hi densities 
measured for this system. This can easily elevate 
Nhi by a factor of 3 (0.5 dex) alone, so the actual 
fraction of gas above the iMavbhate et al.l ( 2007 ) 
threshold (and at high mechanical energies) may 
be small. Consequently, the idea of a critical Hi 
column density for massive cluster formation may 
still be supported in this case. 

A.2. NGC 1487W 

This tail has no SCCs according to the present 
optical criteria. Like the eastern tail, optical de- 
tection of point sources in Mil demonstrate that 
this tail may preferentially exhibit a star form- 
ing morphology devoid of massive clusters. In- 
stead, this tail hosts a large number of fainter 
blue (y - / - 0.0) sources across the WFPC2 
FOV (Mil). The distributions of Nri, ctios, and 
Eke are very similar to those of this tail's eastern 
counterpart, though again it is likely that projec- 



tion effects may artificially enhance the apparent 
column densities and affect other measured quan- 
tities. Overall, dvios/dr± in this tail appears con- 
strained to values generally <20 km s^^ kpc~^, 
unlike in NGC 1487E, which shows several regions 
with larger velocity gradients. Otherwise, ^kpc- 
scale mechanical energies across the Hi debris are 
similar between these two tails. 

A.3. NGC 4747 

The debris seen in NGC 4747 results from 
a sa320 Myr old merger between its pr ogenitor 



galaxy and the more ma. ssive NGC 4725 ( Havnea 
19791 : IWeversetal.lll984l) . As in NGC 1487E/W, 
sources were observed at fainter magnitudes in 
Mil, but few sources (6 in tail) were detected at 
My < -8.5. These also have a good probability 
of being real clusters based on star counts; /sc = 
0.80 overall, while one pixel has a Preai value >0.94 
evident in Figure 5.3. 

The WSRT data used in this analysis is not 
as deep as the VLA/ATCA cubes for other tails, 
but accounting for the beam dimensions, it is still 
sensitive with our channel masking to log Njji 
> 19.6 (cm~^). A small amount of low column 
density gas is detected towards the center of the 
WFPC2 FOV, and has relatively low velocity dis- 
persions (<30 km~^) on its «6 kpc beam scale. 
We measure high velocity gradients across this Hi 
gas {dvios/dr± > 30 km s~-^ kpc~-^), but the few 
surrounding pixels for each measurement makes 
this quantity prone to uncertainty. A large frac- 
tion of pixels (~50%) are marked as unfitted for at 
least one of the Hi properties plotted, appearing 
in optical images but showing no clear Hi emis- 
sion here. This tail appears optically prominent 
but relatively Hi-poor. 

The fact that this intermediate-stage merger 
has hardly any Hi to the best of its cube's offered 
resolution is peculiar. It may be possible that it 
has been stripped out and/or dispersed by its re- 
cent encounter with NGC 4725. What gas remains 
is characterized by low kinetic energy densities 
(~10*^ erg pc^^) and low column densities. How- 
ever, the compliment of blue {V — I < 0.5) point 
sources seen in Mil may indicate that less mas- 
sive star clusters (^10^ -^0) have formed in the 
past few hundreds of Myr. Thus, it appears that 
the formation of these structures would require a 
lower column density threshold than the fiducial 
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log Nhi fti 20.6 (cm ^) proposed for more massive 



clusters (jMavbhate et al.M2007r ). and/or the star 



formation that has occurred may have helped effi- 
ciently consume/disperse an existing Hi reservoir 
over the past few hundred Myr. 

A.4. NGC 520 

The tip of NGC 520 tail, the Hi and op- 
tical condensation seen here and in Mil could 
be a TDG from a wSOO Myr old interac- 
tion or a d warf galaxy entrained i n the Hi 
tidal stream (jStanford fc BalcelldlliQlh . Within, 
fainter sources than the SCC magnitude cutoff 
were originally detected. Two sufficiently bright 
objects of the Mil source lists an be considered 
SCCs and are scattered in the Hi debris, but they 
have low probabilities of being real clusters (fsc 
= 0.06). 

The very center of the potential TDG contains 
high column density gas (Nhi > lO^'^'^ cm^^) on 
w3 kpc scales, and has a low linc-of-sight velocity 
gradient {dv\os/dr± < 5 km s^^ kpc^^). However, 
the Ske distribution peaks barely above 10^^ erg 
pc~^ from the moderately turbulent nature of the 
gas, and overall the fraction of the tail Hi that is 
both of high column density and turbulent on sa3 
kpc scales is ^17%. Once more, this environment 
may be somehow favorable for a smoother, diffuse 
star formation morphology that does not permit 
the formation of massive star clusters. 

A.5. NGC 2992 

NGC 29 92 is ftilOO Myr in to its interaction with 
NGC 2993 (|Duc et allbOOOh . It is heavily inchned 
along our line of sight, enhancing the tail column 
density by a factor of ~3 (0.5 dex). The putative 
dwarf galaxy candidate at the tail tip was previ- 
ously confirmed to be kinema tically distinct and 
therefore is likely a re al TDG (JBrinks et al.ll2004 : 
iBournaud et~aLl |2004| ) . Most of the 10 SCCs de- 
tected in this system can be found there, at least 
in projection. Though Mil found that across the 
debris, NGC 2992 is deficient in clustered star for- 
mation (to My < -8.5), individually we find that 
about 72% of the detected sources are probably 
massive clusters, especially given their nonuniform 
distribution. 

The Nhi measurements for in-tail pixels are all 
high in projection, with almost 80% of the gas 



having apparent (w4 kpc-scale) column densities 
log Nhi ^ 20.6 (cm"'^). This is also evident in 
the Smh, rnap, where regions with mass densities 
>10 Mq pc~^ are apparent, despite the satura- 
tion of Hi known to occur around that (dcpro- 
jected) density. We detect a reasonably uniform 
(Tios distribution, with values ranging from Ri20-80 
km s~^, while our map of dvios/dr± highlights a 
region with a locally shallow (<5 km s~^ kpc~^) 
velocity gradient within the TDG and its SCCs. 

In terms of the SCCs, K-S tests reveal that 
the distributions of tail pixels with and without 
at least one SCC are similar for all properties ex- 
cept Nhi, owing to the cluster candidates situated 
in the high column density TDG (though note the 
small sample size for the former category). Again, 
though it is not clear from Figure 5.5, these SCCs 
appear superimposed over a region of especially 
shallow large-scale velocity gradient in Figure 4.5. 
And even if the Nhi values were decreased by 0.5 
dex to roughly account for tail geometry and ori- 
entation. Figure 5.5 would show that these SCCs 
would still be embedded in Hi gas with higher 
column densities than the iMaybhate et al.l (|2007l ) 
threshold. 

A.6. NGC 2993 



The interacting, similarly massed counterpart 
to NGC 2992, NGC 2993 presents a relatively face- 
on orientation. High Nhi areas on «4 kpc scales 
also register high moment-measured velocity dis- 
persions (fJios ^ 40 km s^^) and kinetic energy 
densities (log Ske ^ 47 erg pc~^), and there ap- 
pears to be a strong dependence of aios on Nhi 
in Figure 5.6. About 70% of these sources should 
be clusters, but only half lie in areas of high Nhi- 
Most are still found in regions of high kinetic en- 
ergy density (log Eke ^ 46 erg pc~^). However, 
K-S tests show that the distributions of SCC- 
empty and SCC-laden pixels are not independent 
for any measured quantity. 

A.7. NGC 2782E 

A w200 Myr-old minor merger, NGC 2782 ex- 
hibits an eastern tail that may be the dist orted 
remnant of the smaller galaxy ( Smithj ll99J) . In- 
terestingly, the optical tail extends beyond the Hi 
tail. Approximately 57% of the sa4 kpc-scale tail 
pixels exhibit Hi of supercritical column densities 
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and turbulent velocity dispersions. Many pixels 
are unfitted (/unfit ~ 0.24), so the fraction of high- 
column density, turbulent gas relative to the total 
Hi reservoir is higher. 

The velocity field of the Hi debris varies 
smoothly across the eastern side of the debris 
{dvios/dr± < 5 km s^^ kpc~^), and changes more 
sharply on the western side {dvios/dr± > 10 
km s~^ kpc~^). Velocity dispersions can also 
range from ^10 km s^^ on the eastern side to 
>60 km s^^ in the west. The Ske map follows 
this trend accordingly. We compute a high /gc 
value (0.96) for the 56 SCCs in this tail; thus 
most should be clusters. These nearly all populate 
the turbulent-high Nhi region of Figure 5.7. They 
also lie in pixels with log Ske > 46.5 erg pc~^, and 
span a range of dvios/dr±^ measurements (between 
0-30kms-i kpc^i). 

The pixels with SCCs also comprise different 
distributions of Nhi, ctios, Ske, dvios/dr±, and 
Eke, with K-S tests producing results <5x lO""* 
in all cases. Cluster formation may occur in 
regions of preferentially higher column density, 
higher line-of-sight velocity dispersions, high ki- 
netic energy densities, and a broad range of ^10 
kpc-scale velocity gradients in this tail. These gra- 
dients may be too broadly constructed (measuring 
differences over areas larger than the radio beam) 
to accurately trace the kinematic structure of lo- 
cal SCC-conducive environments, however, partic- 
ularly where the velocity field varies rapidly to- 
wards the central galaxy. 

A.8. NGC 2782W 

If NGC 2782E is the consumed minor galaxy 
of the NGC 2782 interacting pair, this western 
tail may be the Hi-ri ch tidal deb r is of the pri- 
mary galaxy. While ISmith et al.l (|l999( ) detect 



CO emission in the eastern tail, none is found 



in NGC 2782W by .Braine et all (|200l[ ). Per- 



haps correspondingly, far fewer SCCs (6) are de- 
tected in this tail. However, as mentioned in Sec- 



tion mH the work of Torres-Flores et al.l (|2012[ ) 
and our star counts imply that ~52% of these 
sources should be star c lusters. In particular, 
Torres-Flores et al.l ( 20121) find these sources are 
wlO'* Mq, « a few Myr old clust ers. This is con- 



sistent with Ha observations of Knierman et al 



(|2012l ) and their conclusion that ongoing star for- 
mation does occur in this tail, but at low effi- 



ciency compared to other extragalactic environ- 
ments. For no properties do the distributions of 
SCC-carrying and SCC-devoid pixels appear sta- 
tistically distinct. 

The central debris shows large regions with 
very shallow velocity gradients {dvios/dr± < 5 
km s~^ kpc^^), and generally a less extensive 
range of line-of-sight velocity dispersions (cios ~ 
10-30 km s~^) on «4 kpc scales. Compared to 
NGC 2782E, This tail has qualitatively similar gas 
fractions in the turbulent -high Nhi quadrant of 
Figure 5.8. The western tail has a lower fraction 
of unfitted pixels, however, (/unfit ~ 0.1), and a 
higher percentage of gas with relatively low val- 
ues of dvios/dr±. The distributions of Nhi, cios, 
and Eke also peak at broadly lower values for 
NGC 2782W than for the eastern debris. 

A.9. NGC 2444 

A young («100 Myr old) tail, the debris 
of NGC 2444 formed quickly from a pos sible 
high-speed encounter (JAppleton et al.l Il987l : c.f. 
Koopmann et al.ll2008l ). Gas profiles indicate the 
Hi is kinematically disturbed ((Tios~ 20-80 km s~^ 
for most of the pixels) on sa8 kpc scales. There are 
no pixels with Nhi above the fiducial threshold, 
and the dvios/dr± map of Figure 4.9 shows the 
velocity gradients across the tail are mostly below 
wlO km s-^ kpc"^ 

5 of the SCCs detected lie within the tidal de- 
bris; «3 of them may be star clusters based on 
star counts {fsc = 0.63). Moreover, the pix- 
els upon which they appear superimposed do not 
seem to have different distributions of Nhi, cios, 
dvios/dr±, or Eke than do all tail pixels (to the 
extent to small sample of SCC-containing pixels 
allows). Most («95%) of the ^8 kpc-scale Hi has 
kinetic energy densities above 10^^ erg pc~^ and 
has a relatively shallow velocity gradient, but does 
not have a widespread cluster population. If there 
are indeed a few clusters in the SCC source list 
for this system, either they did not need kpc-scale 
high column density Hi gas to form, our resolu- 
tion is too poor to see more local concentrations 
of Hi, or previous star formation dispersed the gas 
efficiently in the past 100 Myr. Given the opti- 
cal faintness of this tail (it docs not register on 
WFPC2 images), it may be unlikely star forma- 
tion feedback has played much of a role in this 
case. 
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A.IO. NGC 2535 

The approximately face-on NGC 2535 debris re- 
sults from a younR (ftilOO Myr) interaction with 



NGC 2536 (jHancocket al.ll2007n . Measurements 
of «6 kpc-scale velocity dispersions for these pro- 
files ranges from 10-50 km s~^, typically, with the 
regions of highest aios and lowest dvios/dr± (< 
4 km s^^ kpc^^) closely mirroring the most opti- 
cally prominent regions seen in Mil. Most of the 
36 in-tail SCCs also fall in these regions, «93% 
of which should be star clusters. From K-S tests, 
the pixels that contain them appear to have gen- 
erally higher values of Nhi, and lower values of 
dv\os/dr^, but do not arise from different distri- 
butions of (Tios and Ske compared to all tail pix- 
els. Several pixels have multiple SCCs "within," so 
their p^-caX values are higher (>0.94); these occupy 
the regions of turbulent, high Nhi, low dv\os/dr±, 
and high Ske in Figure 5.10. 



A.ll. NGC 6872E 

NGC 6872 is sal50 Myr into an 
action with t he smaller companion 
IC 4970 (iHorellou fc Koribalski 



inter- 
alaxy 



20031 



Horellou fc Koribalskil |2007|) . NGC 6872E is 
the easternmost FOV of the enormous eastern 
tail, featuring a moderate percentage of high 
column density Hi (about 50% of the gas has 
log Nhi > 20.6 cm^^ on «16 kpc scales). The 
line-of-sight velocity dispersion ranges from «20 
km s~^ from the tail edges to «100 km s~^ in the 
interior. 

The central region of the debris is also host 
to an extensive population of SCCs (over a hun- 
dred are detected here), which is also coincident 
with areas of shallow large-scale velocity gradients 
{dvios/dr± < 5 km s~^ kpc"'^). For all considered 
properties except dvios/dr±, the pixels hosting 
SCCs comprise different distributions than those 
of tail pixels that do not; their distributions tend 
to peak at higher values of Nhi, cios, and there- 
fore also Eke- According to Table 3, about 94% of 
these sources should be clusters, and pixels with 
multiple SCCs have even higher probabilities of 
containing at least one. As seen in Figure 5.11, of 
these pixels of highest cluster-containing probabil- 
ity, almost all have kinetic energy densities >10^^ 
erg pc~^, and column densities >10^°'^ cm~^. The 
positions of the SCCs are also aligned with regions 



of Ha emission ( Bastian et al.ll2005t iMihos et al 
1993). 

There are still cluster candidates in pixels with 
"subcritical" Nhi, however — including roughly 6- 
10% of them that lie in regions unfitted by several 
Hi properties — and many of these may be real. 
It is then again a question of whether this Nhi 
threshold is important, these clusters had formed 
in previously high column density Hi in the past 
150 Myr that has since dispersed or been con- 
sumed, or instrumental effects prevent an assess- 
ment of Hi properties at crucially lower (pc-kpc) 
scales. A non-negligible fraction of the tail pix- 
els are unfitted (/unfit ~ 0.15), so these effects are 
certainly possible. For NGC 6872, the physical 
scale may be especially important, as the cube for 
this system has a particularly large beam, both in 
angular and physical extent. 

A.12. NGC 6872W 

This western debris highlights the base of the 
eastern tail of the NGC 6872 system, displaying a 
central clump of moderate Hi emission around the 
Mavbhate et al.l ( 20071 ) limit, immersed in a tail of 
lower column density gas (~80% of the in-tail pix- 
els) that follows the optical structure seen in Mil. 
Flanking the Hi tail is bright optical debris that 
both contains SCCs and is devoid of detectable 
neutral hydrogen. This helps account for the large 
fractions of unfitted SCC-hosting and SCC-empty 
pixels («0.4 and 0.3, respectively). The densest 
regions exhibit velocity dispersions up to «100 
km s~^. Here the velocity structure is not of a 
smoothly rotating, extended disk, but has a num- 
ber of kinematically distinc t components, as elab- 
orated in the modeling of iHorellou fc Koribalski 
POOTI ). Thus, the dvios/dr± map, showing mea- 
surements of line-of-sight velocity gradients across 
pixels separated by more than one beam, displays 
areas with seemingly steep gradients. Internally, 
however, the local gradients may be shallower. 

In either case, the 101 SCCs detected appear 
as widely distributed as the kinematically varied 
Hi debris. K-S tests demonstrate that the pix- 
els that contain them are generally not different 
in terms of Nhi, ctios, or Eke, and while their dis- 
tribution of dv\os/dr± values appears statistically 
distinct, the aforementioned spatial resolution of 
these NGC 6872 data make this difficult to inter- 
pret. 18 of these sources exist in exclusively op- 
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tically defined debris, and the remaining sources 
mostly in Iow-Nhi areas are also likely to be actual 
clusters (/sc = 0.93). Even the pixels with the 
highest probabilities of having at least one cluster 
generally have Nhi values below the fiducial cut- 
off, while their high ctios measurements ensure that 
their mechanical energy densities remain all above 
log Ske = 46 ergpc"^ 

A.13. NGC 1614N 

The northern region of the NGC 1614 tidal de- 
bris is Hi-poor (/unfit ~ 0.89), with concrete de- 
tections only at the very northern and southern 
edges of the central host galaxy. It is conversely 
prominent in the optical, with galactic material 
arranged in warped rings around the merger rem- 
nant of a n interaction thou ght to have begun «750 
Myr ago ( Neff et al.lll990l and references therein). 
The few pixels with detectable Hi emission show 
high (Tios (~75 km s~^). These isolated clumps 
of neutral hydrogen have high column densities 
on «5 kpc scales, and appear to have larger-scale 
velocity gradients dv\os/dr± '-^lO km s~^ kpc""'^, 
though the paucity of Hi pixels makes this com- 
parative measurement admittedly uncertain. 

19 SCCs are found interspersed in the optical 
debris, with only 3 lodged in the Hi debris (sev- 
eral occupy the same pixels). Thus, K-S tests can- 
not be performed on the different pixel distribu- 
tions. But the majority of the SCCs should be 
clusters {fsc— 0.95), so clusters are detected in 
this tail regardless of Hi content. NGC 1614N re- 
sults from a very late-stage interaction, so these 
sources may also be several hundred Myr old and 
had initially formed while surrounded by an Hi 
medium that has since dispersed, been ionized, or 
consumed. While the photometry of Mil cannot 
age-date the SCCs, their red V — I colors {V — I^ 
1) may imply these clusters are older (or reddened 
by extinction). 

NGC 1614 is luminous at infrared wave- 
lengths and conta ins a powerful nuclear starburst 
(JAlonso-Herrero et al. 2001). It is possible that 
the Hi in the northern tail has been eliminated by 
central winds, photoionized by the central star- 
b urst, or pho t oioniz ed by past star formation. 



Hibbard et al. 



(2000|) contends that a mean R- 
band surface brightness of «23 mag arcsec"^ — 
depending on the Hi column density and dimen- 
sions of the tidal debris — is required for the stellar 



population to effectively rid itself of the warm neu- 
tral medium on these scales. NGC 1614N is one of 
the few debris fields in the Mil sample to exceed 
this surface brightness (its F606W-band surface 
brightness is «22 mag arcsec"^, which in terms 
of photometric bandpasses is close to R), so it is 
possible mass in situ photoionization of the tidal 
Hi has occurred in the past. 

A.14. NGC 1614S 

The Hi of the southern tail of NGC 1614 is 
confined to a "stub" extending from the southern- 
most tip of the merger remnant; most tail pixels 
are optically defined (/unfit ~ 0.89). Its Hi char- 
acteristics are similar to those of the northern de- 
bris; the measurable Hi here is of comparably high 
column density (log Nhi ^ 20.6 cm~^) and with 
line-of-sight velocity dispersions ranging from the 
instrumental limit to ~75 km s^^. The ranges of 
velocity gradients are also similar, with dvios/dr± 
« 5-15 km s~^ kpc~^. 

The main difference is the number and distri- 
bution of SCCs. 30 sources are located in the com- 
bined optical and Hi tail, with most of them falling 
in the Hi stub. Once more, a high fraction of them 
are probable star clusters (/sc = 0.96, with psc 
exceeding this number for many pixels that con- 
tain more than one SCC). The sources here are 
very blue {V — I < 0; Mil), indicating star for- 
mation might be relatively recent and/or ongoing. 
These sources trace out the high-Nni, high-SKE, 
and turbulent regions of the Hi parameter space of 
Figure 5.14, and the pixels that have them do not 
appear different with respect to those properties 
than Hi pixels overall. 

It is however curious that the positions of these 
SCCs is offset by a few arcseconds from the re- 
gions of highest column density as seen in the Nhi 
map of Figure 4.14. The Hi may be in the process 
of being locally consumed or swept out by current 
star formation. Approximately 42% of the equally 
high probability SCCs are found in optically bright 
yet Hi-unfitted regions, so these sources may be 
further examples of past feedback and H i disper- 
sal or consumption sometime over the course of 
this tail's «750 Myr lifetime. Unfortunately, VI 
observations of Mil do not permit accurate age- 
dating of cluster candidates to provide evidence for 
this age-dependent Hi removal. If generous K-S 
probability cutoff is employed, the distributions of 
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dvios/dr± between pixels and without SCCs ap- 
pear marginally independent. But again the low 
spatial resolution for this especially distant object 
may confuse this quantification of the large-scale 
velocity field. 

A.15. NGC 4038A/B/C 

This well-studied merger has two large, sym- 
metric tails of ages r:!450 Myr (K03 and references 
therein). WFPC2 images of the ABC positions 
along these tails reveal them to be optically faint, 
diffuse, and all but devoid of SCCs (K03; Mil). 
As a whole they are poor in Hi, w i th ver y little gas 
at or above the iMavbhate et al.l ( 20071) threshold 
on wkpc scales. Across the WFPC2 pointings, it 
appears «30-80% of the tail pixels are classified 
as unfitted by at least one Hi characteristic. This 
gas is relatively kinematically unperturbed, with 
cios ^ 25 km s~^, and much of the gas showing 
dispersions less than the turbulent-thermal cutoff. 
Consequently, less than r:!5-20% of the gas has ki- 
netic energy densities >10*^ erg pc^^. 

3 SCCs are found in NGC 4038A, and 1 is found 
in the NGC 4038B and NGC 4038C FOVs each. 
They all have low values for /sc and their host 
pixels have correspondingly low values of Proai (0- 
0.24); thus they are are likely point source con- 
taminants. They populate too few pixels to con- 
duct K-S tests between distributions of Hi char- 
acteristics. Unlike the active, luminous interior 
and bridge regions of the merger, the tails of this 
system are extraordinarily quiescent in terms of 
ongoing star formation. 

A.16. NGC 3256E/W 

Dynamically, the tails of NGC 3256 are about 
400 Myr old (K03), resulting from an ongoing ma- 
jor merger of two similarly massed galaxies. The 
channel width of the NGC 3256 data cube is wide 
(«33.5 km s""'^), so each pixel has few channels of 
signal. Combined with the beam size and sensi- 
tivity of the data, uncertainties in Hi properties 
are high for pixels with Nhi S 10^° cm~^. The 
vertical "tail" of constant a\os in Figures 5.18 and 
5.19 represents the instrumental limit in determin- 
ing velocity dispersions. Much of the Hi gas in 
these tails is of high column density and line-of- 
sight velocity dispersions on scales of ~5 kpc; we 
find «72% and 88% of tail pixels has those proper- 



ties in the eastern and western tails, respectively. 
The distribution of dvios/d,r± values in both tails 
is very broad, ranging from 0-50 km s~^ kpc~^, 
while the distribution of that property is skewed 
towards lower values for NGC 3256W. 

The western tail has fewer pixels with at least 
one unfitted Hi characteristic (by a factor of sa4- 
10). 32 SCCs are located in the eastern tail, and 
51 are in the west. Given the calculated values 
for /sc (0.50-0.60), a little over a dozen sources 
in the east and about 30 sources in the west may 
be genuine star clusters. This strong level of con- 
tamination results from this system's low galactic 
latitude (5 « 11°), and led K03 to conclude that 
the eastern tail did not have a statistically signifi- 
cant widespread cluster population, and the west- 
ern tail did. This is still true; the eastern tail may 
have only a handful of real clusters, but they are 
spread throughout the sal 720 kpc^ tail area and 
do not constitute a significant population. 

Most of the SCCs in either tail populate the 
high Nhi and high crios parameter space of the 
neutral medium. According to K-S tests, pixels 
containing SSCs in the eastern tail may have typ- 
ically higher line-of-sight velocity dispersions, but 
do not comprise unique distributions compared to 
all pixels in terms of other measured properties. 
Conversely, the pixels hosting cluster candidates 
in NGC 3256W might have generally shallower 
velocity gradients with a generous interpretation 
of that K-S result. The central "spine" of the 
western debris is a region of low dvios/dr± (<5 
km s~^ kpc"-*^), an attribute absent in the eastern 
tail. Perhaps a combination of this and the higher 
fraction of high column density gas contributes to 
the factor of ^3 increase of possible clusters in the 
western tail over the east. 

A.17. NGC 7252E/W 

NGC 7252 is a late-stage (w730 Myr old; K03) 
merger between two comparably massed galaxies. 
The channel width of this system's data cube is 
the widest of the debris sample (w42.5 km s~^). 
The Hi in both debris fields is relegated to the 
tips of the straight and narrow eastern and west- 
ern tails, and is characterized by low column and 
mass densities on « 1.5 kpc scales and steep veloc- 
ity gradients {dvios/dr± > 15 km s~^ kpc~^) on 
larger scales. The latter may be attributed in part 
to the velocities ascribed to the low-signal edge 
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pixels, which can be uncertain. The moderately 
turbulent kinematics of the gas helps to offset the 
low mass densities and maintains the kinetic en- 
ergy densities of much of the Hi above 10^^ erg 
pc . 

Comparable numbers of SCCs are in either tail 
(11 and 17 for east and west, respectively), with 
fairly comparable values of /sc (0.80 and 0.89). In 
7252E these sources are widely distributed, many 
of which are situated in the Hi-poor optical debris 
(/unfit ~ 0.40). Two neighboring sources are lo- 
cated in a clump of diffuse optical emission, which 
is slightly offset from the central high-Nni region 
seen in Figure 4.20. Pixels that contain the SCCs 
in this tail do not seem to arise from dissimilar 
distributions of Nhi, tios, Ske, or dvios/dr± than 
all Hi tail pixels. 

The same is true for the western tail, except for 
the velocity dispersion, for which SCC-laden pixels 
seem to have a more narrow range of values (ctios ~ 
10-30 km s~^). 11 sources are dispersed through- 
out the primarily optically defined tidal tail, with 
«40% in especially Hipoor regions (in comparison, 
/unfit ~ 0.30). A few Hi pixels in the only high 
mass density region at the end of the tail contain 
6 of them. This region of the debris has the lowest 
dvios/dr±, and is the only part of Figure 5.21 that 
shows Nhi > 10^°'^ cm^^. These sources belong 
to the promi nent optical/Hi-defin ed TDG in the 



western tail ( Hibbard et al.lll994 ). The fact that 



this one region of obvious clustered star formation 
is the only region of high column density Hi and 
contains a large number of SCCs (the superim- 
posed black points in Fig ure 5.21) indicates that , 



at least in this case, the iMavbhate et al.l ([20071) 



Nhi threshold may be important in forming star 
clusters. These few Hi pixels do not stand out 
in the distribution of SkeJ & combination of high 
mechanical energy densities and column densities 
may be observationally important in finding SCCs 
in situations like this. 



completeness limit at My = -8.8. Therefore, it 
is probable that some SCCs are missing from the 
source list. 

Its Hi debris has velocity dispersions on this 
scale of ~5-30 km s~^, and has a smoothly vary- 
ing velocity structure with shallow larger-scale ve- 
locity gradients {dv\as/dr± < 5 km s^^ kpc^^). 
Approximately 10% of the tail pixels has a high 
column density, and is confined to the long tail- 
shaped structure that extends from the north to 
the south of the Nhi map. Of the 14 SCCs in 
the Hi or optical debris, only two are embedded 
in this relatively higher- Nhi material, and most 
are dispersed throughout the debris. About 87% 
of these sources should be clusters, and one pixel 
that contains two sources — part of a short chain 
of sources that extends from the north of the op- 
tical image — correspondingly has a high chance of 
containing at least one massive cluster. 

These SCC-supporting pixels do not appear to 
have statistically unique distributions of measured 
Hi quantities compared to all Hi pixels, and all 
but one of them lie in regions of low column den- 
sity (log Nhi < 20.6 cm~^). Their moderate ve- 
locity dispersions mean that their kinetic energies 
are still greater than the fiducial 10^^ erg pc~^. 
As is the case with other tails, the low spatial res- 
olution, feedback from star formation (certainly 
possible over the w460 Myr lifetime of the tail; 
K03), and/or a revised Nhi threshold might at 
least partially explain the detection of these hand- 
ful of reasonably probable clusters. The particu- 
larly poor resolution makes it especially possible 
that sub-kpc clumps of denser Hi gas are being 
effectively diluted by the size of the beam for this 
system, especially when they are found in closer, 
better resolved tails. 



A.18. NGC 3921S 

NGC 3921 is the most distant interaction of the 
sample (84.5 Mpc), and accordingly its Hi cube 
has the worst spatial resolution of the debris sur- 
veyed here («2 kpc for each pixel; roughly 4 times 
that for the size of the beam). The optical com- 
pleteness limit for this tail as analyzed by K03 is 
also the worst of the sample, with a y-band 50% 
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